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by 
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Abstract 

In  order  to  investigate  beam  deformation  and  depletion  at  very  high 
photon  fluxes,  a  computational  procedure  has  been  developed  that  is  suitable 
to  simulate  the  propagation  of  intense  beams  through  transparent  solids. 

It  was  applied  to  NaCl  exposed -to  a  FWHM  21  psec  pulse  from  a  frequency- 
doubled  Nd-YAG  laser  (A  =  0.532  ym) .  The  energy  content  of  this  TEMqo~ 
mode  pulse  is  equal  to  the  experimentally  determined  single-shot  damage 
threshold  of  this  material.  The  most  convenient  mode  of  data  display  was 
found  to  be  a  motion  picture  showing  the  propagation  of  the  pulse,  represented 
as  a  colored  contour  map  of  the  local  photon  flux,  through  the  focal  region. 
Similar  displays  of  the  spatio-temporal  development  of  the  free  carrier  con¬ 
centration  have  been  made.  The  material-laser  photon  interactions  in  these 
calculations  were  based  on  the  two  main  laser  damage  mechanisms:  avalanche 
ionization  and  multiphoton-polaron  absorption.  The  most  significant  result 
of  these  investigations  is  that,  according  to  both  models,  beam  deformation 
prevents  the  build-up  of  a  sufficient  photon  flux  in  the  interaction  area 
to  cause  damage  at  this  pulse  energy.  This  means  that  either  both  damage 
mechanisms  do  not  explain  the  best  experimentally  available  data  or  these 
data  are  extrinsic  rather  than  intrinsic  features  of  the  material. 

Spherical  aberration  of  the  focusing  lens  as  a  possible  cause  for  the 
laser  damage  morphology  observed  at  threshold  in  NaCl  has  been  discussed 
as  well.  Independently,  the  response  of  single-crystalline  germanium 


exposed  to  intense  pulses  of  0.45  eV  photons  was  theoretically  studied. 
Quantitative  calculations  of  the  relevant  multi-photon  cross  section  for 
free  carrier  generation  and  of  the  dielectric  response  function  have  been 
performed.  These  results  have  been  coupled  with  a  model,  taken  from  the 
literature,  which  describes  the  resulting  changes  in  the  energy  distribution 
of  the  electrons,  holes  and  both  acoustical  and  optical  phonons.  These 
studies  represent  the  first  detailed  theoretical  calculation  of  the  laser 
damage  threshold  of  germanium  at  this  photon  energy. 
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I.  INTRODUCTION 


Exposure  of  transparent  materials  to  intense  laser  beams  induces 

changes  of  their  optical  properties  (conveniently  expressed  in  terms  of  the 

dielectric  function  s(u)).  This  in  turn  influences  the  photon  trajectories 

of  any  beam  whose  intensity  profile  is  not  uniform.  The  mutual  interaction 

of  the  photon  field  and  the  solid  is  very  complex.  The  majority  of  attempts 

to  calculate  the  spatio-temporal  field  distribution  in  a  transparent  solid 

having  a  field-dependent  dielectric  function  are  essentially  based  on  solutions  of 

Maxwell's  equations. In  these  calculations  it  is  commonly  assumed  that 

2 

the  linear  refractive  index  simply  increases  by  ^E  in  the  presence  of  an 

optical  field  of  rras  field  strength  E.  Values  for  the  so-called 

nonlinear  refractive  index  ^  are  taken  from  measurements  (e.g.  Smith^) 

or  independent  calculations.  The  most  pronounced  effect  produced 
2 

by  the  term  n^E  is  self-focusing  of  laser  beams  having,  e.g.,  a  Gaussian 

intensity  profile  (TEM  mode).  A  characteristic  feature  of  these  field 

oo 

theories  of  beam  deformation  is  a  singularity  of  the  field  strength  in  the 
focal  waist.  As  a  result,  the  distribution  of  the  rms  optical  field  strength 
cannot  be  determined  in  the  space  behind  the  focal  plane.  Attempts  to  do 
so  invariably  violate  the  law  of  conservation  of  energy.  ^ 


Little  if  any  work  has  been  done  on  beam  deformation  due  to  the 

reduction  of  the  refractive  index  caused  by  the  production  of  photo  carriers 

(u  5)  (*^  d  ^ 

in  the  intense  photon  lieid.  ’  Mulciphoton  absorption  or  electron 

(2) 

impact  ionization  of  lattice  constituents  are  mechanisms  for  the  genera¬ 
tion  of  free  carriers.  Experimental^'^  as  well  as  theoretical^  evidence 

21  -3 

exists  for  free  carrier  concentrations  to  reach  up  to  10  cm  at  pnoton 
fluxes  corresponding  to  tne  dielectric  breakdown  threshold  of  the  solid. 


It  is  likely  that  this  negative  contribution  to  the  refractive  index 
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2  17  -3 

to  dominate  n,E  already  at  concentrations  around  10  cm  or  less  in 

alkali  halides  exposed  to  short  laser  pulses  of  optical  frequencies.  '  ’ 

This  clearly  indicates  that  beam  deformation  (self-defocusing)  via  a 

reduction  of  the  refractive  index  is  quite  possibly  occurring  at  power 

densities  well  below  the  laser  damage  threshold.  A  better  understanding 

of  this  phenomenon  is  clearly  called  for  in  view  of  the  fact  that,  during 

the  operation  of  high  power  lasers,  power  densities  approaching  and — nil 

too  frequently — exceeding  the  damage  theshold  of  exposed  optical  materials 

must  routinely  be  handled.  In  addition  to  these  largely  economical 

reasons,  unknown  self-focusing  and  defocusing  effects  prevent,  for  example, 

the  exact  measurement  of  the  damage  threshold  of  a  transparent  solid  that 

Is  subjected  to  a  short  laser  pulse.  Such  experiments  are  commonly 

executed  by  tightly  focusing  the  laser  beam  into  the  bulk  cf  the  sample. 

Since  it  is  impossible  to  directly  measure  t  •  >  photon  flux  In  the  focal 

region,  it  is  calculated  on  the  basis  of  diffraction-limited  optics,  corrected 

(2) 

in  first  order  for  self-focusing  only.  Photon-induced  beam  deforma¬ 
tion  caused,  e.g.,  by  negative  contributions  to  the  refractive  index,  has 

(9) 

so  far  never  been  considered  in  any  detail  in  the  context  of  laser  damage. 

An  exception  is  an  interesting  piece  of  information  on  the  importance  of 
beam  deformation  that  was  provided  by  Danileiko,  Manenkov  and  Sidorin^^ 
who  tried  to  damage  Gs  at  A  =  2.76  urn  and  A  =  2.94  ;;n  with  nsec  pulses. 

These  workers  did  not  succeed  because  oL,  ns  they  speculate,  f roe-carrier 
induced  beam  sel f -defocusing.  No  specifics  were  given. 

A  further  aspect  of  beam  deformation  are  spherical  aberrations  caused 
by  the  focusing  optics  used  in  typical  laser  damage  experiments.  The  peculiar 
damage  morphology  observed  at  damage  threshold ^ 1 1 ^  was  suspected  to  be  in 
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(12) 

part  due  to  this  effect.  Preliminary  calculations  by  Daniieiko  et  al. 
indicated  that  strong  spatial  intensity  oscillations  can  occur  if  lenses 
with  residual  spherical  aberration  are  employed.  We  have  performed  a 
detailed  study  of  these  effects  (see  Chapter  IV)  and  concluded  that  the 
spacing  between  the  obtained  intensity  peaks  does  not  correspond  to  the 
observed  spacing  between  damage  sites  in  the  focal  volume  along  the  laser 
beam  axis.  Yet,  a  considerable  reduction  in  the  peak  photon  flux  may  occur 
in  lenses  with  slight  aberration  and,  thus,  influence  interpretation  of  the 
damage  experiments.  It  was  the  goal  of  the  work  performed  under  this  con¬ 
tract  to  contribute  toward  our  knowledge  of  the  interaction  of  intense 
laser  beams  with  nominally  transparent  optical  materials.  The  approach 
chosen  to  achieve  this  is  outlined  in  the  following  section. 


f 
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II.  RESEARCH  OBJECTIVES  AND  APPROACH 

The  purpose  of  the  work  carried  out  under  this  contract  was  to 
theoretically  study  several  aspects  of  the  interactions  of  very  intense 
laser  pulses  with  solids.  Emphasis  was  to  be  placed  on  infrared  photons 
and  infrared  optical  materials. 

In  particular,  we  have  investigated  the  changes  of  the  dielectric 
response  function  due  to  these  interactions  and  the  resulting  deformation 
and  attenuation  of  the  laser  beam.  Aberration  effects  induced  by  spherical 
focusing  lenses  have  been  discussed  as  well. 

These  goals  required  to  not  only  consider  the  temporal  (dynamic) 
characteristics  of  the  photon-solid  interactions  but  also  to  include  their 
spatial  variations.  Spherical  aberration  was  studied  within  the  paraxial 
approximation  of  the  wave  equation,  excluding  any  non-linear  interaction 
effects  such  as  self-focusing,  in  an  attempt  to  interpret  the  peculiar 
damage  morphology  observed  in  the  focal  volume  of  the  laser  beam  at  threshold 
peak  flux.  The  nonlinear  aspects  of  the  photon-induced  changes  of  the 
dielectric  function  required  the  adoption  and  development  of  a  special  com¬ 
puter  code  "PULSE".  It  became  operational  in  the  first  contract  period 
(Aug.  1,  1978  -  July  31,  1979)  and  was  first  tested  and  refined  by  simu¬ 
lating  the  classic  experiments  by  Smith  et  al.^^  who  focused  a  21  psec 
FWHM  pulse  at  0.532  pm  into  NaCl  and  determined  the  damage  threshold.  A 
motion  picture  of  this  damaging  pulse,  traversing  the  focal  region  of  a 
2.54  cm  focal  length  lens  inside  NaCl,  was  deposited  at  AFOSR.  Since  then 
we  have  included  beam  depletion  in  the  code  and  refined  the  attainable 
spatial  resolution.  We  applied  it  to  a  comparison  of  beam  deformation  as 

obtained  under  the  above  mentioned  conditions  within  the  avalanche  ioni- 
(  2  ) 

zat ion  mechanism  and  the  mult iphoton-po Laron  ahsoption  model  ot  laser 


damage.  The  results  are  discussed  in  Chapter  III. 
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The  developed  computer  code  is  based  on  a  particle-in  all  method 

(13) 

originally  used  by  Dudder  and  Henderson.  It  is  designed  to  simulate 

the  complete  spatio-temporal  behavior  of  a  short,  intense  pulse  that  is 
focused  into  an  optical  material,  independent  of  the  physical  mechnisms 
responsible  for  the  photon-induced  charges  in  the  refractive  index  and  the 
temperature  of  the  exposed  solid.  Thus,  the  code  is  sufficiently  flexible 
to  accommodate  any  flux,  wavelength,  material  and  temperature  dependence 
of  the  dielectric  function,  any  possible  mechanism  of  electronic  excitation 
of  the  solid  (such  as  multi-photon  photo-carrier  generation,  avalanche 
ionization  and  the  like)  as  well  as  various  conceivable  mechanisms  for 
converting  energy  of  the  photon  field  to  excited  electron  and  phonon  dis¬ 
tributions  and  ultimately,  to  an  increase  in  the  lattice  temperature.  A 
description  of  "PULSE"  is  given  in  Attachment  A.  The  procedure  used  to 
incorporate  beam  depletion  is  discussed  in  Chapter  III. 

As  stated  above  the  calculations  of  beam  deformation  at  the  wave¬ 
length  in  the  visible  were  performed  mainly  to  test  the  computer  code  and 
to  demonstrate  the  feasibility  of  the  chosen  computational  approach.  Diffi¬ 
culties  associated  with  singularities  in  the  focal  waist  or  with  energy  con¬ 
servation,  characteristic  for  the  theories  based  on  solving  field  equations, 
clearly  do  not  occur  in  our  case.  We  have,  for  the  first  time,  simulated 
the  rather  dramatic  break-up  of  the  laser  beam  expected  under  the  chosen 
high-power  conditions  ( the  energy  content  of  the  pulse  simulated  corresponds 
to  the  so-called  single-shot  damage  threshold  which  has  been  shown  to  be 
well-defined  in  NaCl^^). 

In  order  to  apply  the  developed  computational  method  to  materials 
and  lasers  of  interest  for  use  in  the  infrared  region  of  the  electromagnetic 
spectrum,  it  was  necessary  to  establish  the  mechanism  of  photon-induced 


material  modifications  for  this  case.  Quantitative  information  on  cross 
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sections  for  free  carrier  generation,  photon  absorption  by  free  carriers 

and,  most  importantly,  for  the  calculation  of  beam  deformation,  on  the 

photon- induced  changes  of  the  dielectric  response  function  is  required. 

We  have  approached  this  problem  by  initially  selecting  single-crystalline 

germanium  as  a  typical  infrared  optical  material  that  is  exposed  to  an 

intense  beam  of  0.45  eV  photons.  Calculations  of  the  involved  multiphoton 

(2) 

absorption  cross  section  a  (two-photon  absorption),  as  well  as  of  the 
dielectric  function  have  been  performed  for  this  case.  They  are  discussed 
in  Chapters  V  and  VI  of  this  report. 
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REFORMATION  OF  INTENSE  LASER  BEAMS  TIGrTU  'CC’JSEj  INS  IRE  NaCl  : 

A  COMPARISON  OF  THE  HULTIPHOTON-PQURON  AND  AVALANCHE  MODELS  OF  OPTICAL  BPEAnRuHN 


P.  Kelly  and  D.  Ritcnie 
National  Resea rcn  Council 
Ottawa,  Canada  <1  AOS  1 

and 

P.  Braunlich,  A.  Schmid,  and  G.  W.  Bryant 
Department  of  Physics 
Washington  State  University 
Pullman,  Washington  09164 


Computer  simulations  of  the  interaction  of  intense  picosecond  laser  pulses  ('  =  532  nm) 
with  NaCl  in  tne  focal  volume  of  a  2.54  cm  focal  length  aoerration-free  lens  indicate  tnat 
severe  beam  deformation  may  take  place  at  photon  fluxes  corresponding  to  tne  single-snct 
damage  threshold  measured  under  tnese  experimental  conditions.  Tnis  deformation  is  tne 
result  of  a  chance  in  the  dielectric  function  .vnicn,  in  turn,  is  caused  oy  *ree  earners 
generated  via  multi  photon-assisted  avaianene  ionization  and/or  multipnoton  absorption. 

Strong  sel f-defocusing  prevents  local  onoton  fluxes  from  reaemng  values  reouired  to  raise 
the  lattice  temperature  sufficiently  for  damage  to  occur.  We  conclude  that  either  potn  tne 
avaianene  ionization  as  well  as  the  multi pnoton-polaron  aosorption  mecnanisns  of  laser  preak- 
daiwn  are  invalid  in  tneir  present  form  or  that  the  experimentally  determined  damage  tnresn- 
olrfs  do  not  pertain  to  the  intrinsic  properties  of  NaCl. 

toy  words:  Avalanche  breakdown;  3eam  deformation;  Beam  depletion;  Free  carrier  absorption; 
MuTtlchoton  absorption;  Polaron  absorption. 


1 .  Introduction 
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A  build-up  of  very  large  free  carrier  concentrations  (10  -10  cn  )  during  tne  interaction  of 
intense  laser  pulses  with  wide-gap  optical  materials  is  tnought  to  be  a  prereouisite  for  damage 

occurrence  according  to  all  models  of  laser  breakdown  orooosed  so  far  [1,Z].‘  decent  reviews  of  tnis 
subject  were  given  by  Smith  [3]  and  Srawer  [4],  Two  basic  modes  of  carrier  generation  are  mvo.ea: 
avaianene  ionization  of  valence  electrons  after  multipnoton  generation  of  so-called  starting  elec¬ 
trons  [5,6]  and  multipnoton  ionization  [1,7]. 

The  changes  of  the  dielectric  function 


s(u)  *  £[  *  'Ejj  (!) 

and  the  ensuing  deformation  of  'Gaussian  beams,  associated  with  high  carrier  densities,  nave  been  dis¬ 
cussed  by  Yablonovitcn  and  Eloembergen  [5]  based  on  a  simple  Crude  model  of  conductivity  ana  *'ree 
carrier  absorption,  and  by  Scnnid  et  al.  [2]  wno  retain  tne  Drude  conductivity  term  for  arc 

replace  the  phenomenological  free  earner  absorption  contribution  to  Dy  free-e iectron— acoust  tc- 

pnonon  interaction.  Detailed  calculations  of  the  spatio-temooral  *1ux  distribution  in  :ne  *ccal 
volume  of  an  aberration-free  lens  in  Nad  on  the  basis  of  tne  latter  rnecnamsm  anc  bn  tne  -\d  tipnetbn 
ionization-polaron  aosorption  -’oaei  of  laser  damage  nave  oeen  presented  by  *.elly  et  a'. .  /.  Severe 

beam  sel f-defocusi ng  and  even  beam  oreak-up  was  snown  to  ta*e  place  zje  to  tne  r>ee-crme'-  mjuced 
reduction  in  the  refractive  index.  Seam  depletion,  assumed  to  oe  snail,  «as  neoieccea.  "nese  calcu¬ 
lations  were  performed  in  an  attemot  to  explain  the  damage  moronoiogy  ooserved  by  jrr  tn  and 

coworkers  [9]  as  a  possible  consequence  of  the  comoined  effects  of  sel f-^ocusma  via  n„£~  era  sel*- 

defocusimg.  -iere  n„  is  tne  nonlinear  refractive  index  and  E  tne  rns  field  strenctns  of  the  octicai 
pnoton  field.  c 

"he  cnaracteristic  damage  vestiaes  of  about  f  -in  diameter  and  10  /’•  spacing  aienc  tne  bed’-  axis, 
observed  at  danane  tnresnold  and  described  in  reference  i  9],  could  not  oe  simulated  p>  trese  .'del 
calculations,  quite  possibly  due  to  *nsuf**ic:ert  spatial  resolution.  In  :cnnect'on  *t  ;  :• 

interest  to  note  trat  a  recent  ziscussion  o*  aberration  eTfects  an  me  *iu.«  ::  str-n  .t:  jn  *n  t*v  *xai 
volume  of  snort  focal  lengtn  lenses  CO-/;:  net  account  *or  tne  Ja^aie  ••orpnoioo;/  e’trer  :  ‘  '  tn 

et  al.  ’9;  believe  tnat  local  statistical  variation'-,  j*  tne  dens  it/  oT  start  mo  eiecvc"s  ire  me 
reason  cor  tne  damage  'nicrosites.  Iney  discount  any  effects  tm  free  careers  on  tne  jteiectrm  •  .no¬ 
tion  and  consider  only  self -focusing.  A  calculation  dv  tnese  authors  or  tne  :pat:ai  te-ue^at-re  .’re¬ 
tribution,  i.e.  tne  crystal  volume  inside  an  eilipsoidai  i sotemperat ur  '  surface  *  *  =  1-.74  \ 

is  tne  m^ltina  Point  of  NaCl)  and  its  apparent  coincidence  witn  a  severely  damaged  '•"•non  in  tne  •ocai 
volume  o*  the  laser  beam,  is  taxen  as  one  of  'even!  -.moor  rant  clues  for  tne  validity  jr  rre  a.ji  nc'’? 
model  of  laser  breakdown  [3].  Nevertheless ,  in  lignt  o*  earlier  aiscussions  of  tne  '"ui  tip»  oton- 

ass-.sted  avaianene  model  which  indicated  that  tne  fr*e  carrier  density  nay  reacn  l *?r  rico- 
second  pulses  at  damage  thresnoJ.j  'oj,  a  detailed  computation  of  sel  r-aefocusmo  based  on  me  ava¬ 
lanche  model  was  clearly  Jesiraoie. 

We  nave  applied  the  computer  cod  FULCt,”  developed  by  ►eli/  et  al.  Lol,  to  cm--.  „t otir  am: 
also  recalculated  the  mul  ti  pnoton-polaron  case,  because  recent  improvements  *.n  tne  code  now  accj-nt 
*ully  for  bean  depletion.  Thus,  direct  .omoanson  of  jotn  oarace  mechanisms  is  possible,  ’re 


ig.res  in  'rackets 
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obtained  results  are  remarkable:  neither  cf  the  two  basic  models  of  laser-induced  breax-down  predict 
damage  at  the  experimental ly  determined  threshold  because  or  very  pronounced  sel ‘-defocusmg. 

2.  Models  of  Optical  3reaxaown 

The  details  of  the  avalancne  ionization  model  [3, 4, 6, 9]  and  the  mul tipnoton-polaron  aosorntion 
node  1  j.2,7]  of  laser  breakdown  are  availaole  *rcm  the  literature.  Three  aspects  or  tne  onoton  *rter- 
action  with  optical  materials  are  considered  in  these  model  calculations.  'hev  concern  tne  mechan¬ 
isms  of  carrier  generation  and  their  effects  on  tne  dielectric  function  ana  on  lattice  neating. 


i)  Carrier  generation 

In  NaCl ,  exposed  to  532  nm  pnotons,  the  density,  nc,  of  free  carriers  in  the  conduction  Dana 
changes  according  to 

i 4 )  -4 

n  »  «j.n  ♦  n  r  (2) 


in  the  avalanche  model.  The  second  term  in  eauation  (21  supplies  the  initial  electrons.  However, 
avalanche  ionization  quickly  provides  the  majority  or  carriers  [9].  Its  rate  is  given  by 


log  w.  *  3.58  log  E  ♦  8.62  .  v3) 

Here  E  is  the  rms  electric  field  strength  in' units  of  MV  cnT1  and  is  the  avalanche  ionization  rate 

,i  ’ 

in  sec  Only  multiphoton  carrier  generation  is  retained  in  the  multipnoton-polaron  moael .  All 
symbols  and  parameters  used  above  and  from  nere  on  are  explained  in  table  1. 


ii)  Lattice  heating 

In  the  avalanche  model,  lattice  heating  is  obtained  from  the  Orude-type  conductivity.  'he 
energy  absorbed  per  unit  time  and  volume  is 

VT  *  nce2t  E2/mll  *  ~Z})  .  (4) 


The  detailed  mechanism  by  which  electrons  gain  energy  from  the  photon  field  is  often  referred  to  as 
inverse  Bremsstrahlung. 


According  to  Pokatilov  and  Fomin  [II],  "dressed"  electron  (or  polaron )-phonon  scattering  with 
instantaneous  relaxation  of  the  excited  carriers  pack  to  the  lower  eage  of  the  conduction  oano  is 
taken  as  the  mechanism  py  wnicn  energy  is  gained  by  the  lattice  in  the  multipnoton-polaron  mooei  t2,7j. 
It  has  a  cross  section  ;  and  yields 


copT  »  opFfiu  . 


(5) 


Beam  depletion  is  accounted  for  by  subtracting  the  appropriate  numoer  of  photons  in  ail  processes 

involving  the  aosorotion  of  onotons.  While  this  is  straigntforward  for  processes  involving  :  ana 
g.  the  absorbed  numoer  of  pnotons  per  unit  volume  of  crv* tal  materials  in  tne  avalancne  ionizaf.cn 
P 

process  is  calculated  from  the  energy  gained  by  the  lattice  according  to  eauation  :4,  tpgetner  witn 
the  energy  required  to  produce  an  increase  in  n^  during  a  small  time  interval  wnicn,  cnosen  -or  com¬ 
putational  convenience  and  the  desired  numerical  accuracy,  was  Ct  *  D. 03489  psec. 


Hi)  The  dielectric  function 

Having  independently  accounted  for  absorption  as  described  above,  only  changes  in  the 
refractive  index,  n,  need  be  considered: 


n 


1 


1/2 


rut  . 
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'he  ••eal  and  imaginary  parts  tne  dielectric  function  (ea.  11  are  taken  *:ron  tne  Lruje  node!  or  •>cn 
more  accurate  calculations  of  tne  polaron  conductivity  in  the  presence  or  collisions  «itn  acOuSfc 

phonons  [2],  The  sel  f-focus ing  term  n„£~  is  added  pnenomenologica 1 1  y .  -'or  lacx  or  ••ore  preci  se 

information,  we  nave  simoiy  cnosen  tne  Crude  expressions  for  ail  calculations  based  on  the  avalancne 
model : 

,  x  n2[ I  -  4-e2;2n  ,  mn2  (1  •  -2 '  1 
*  0  CO  J 

and 

'II  3  4"eC'nc/‘n  (9  *  ■ 


The  value  for  the  collision  time  *  is  still  sub’ect  to  much  lebate  .4],  *e  j$ed  the  most  recent  me 
available,  stated  Dy  3r*wer  [4j  to  pe  *  *— j  x  -0*'"  sec,  because  * t  is  also  consistent  witn  *re 
choice  t\w  1  m  rade  in  earlipr  calcu  I  a :  ions  Of  mu  i  1 1  photon  polaron  janaoe  expres¬ 

sion  Dy  Gurevich  et  jl.  is  again  taken  *'or  •  It  is  stated  m  reference  2. 


Jt 


3.  Seam  Ceformat ion 


10 


Beam  deformation  is  a  consequence  of  t^e  nonunirom  beam  profile  (e.q.  -1  udussian;  ana  pnoton- 
inaucea  local  changes  of  the  **er>ac:ive  :maex.  ' re  cc --Outer  core  "r'  ..C£"  was  developed  to  striate 
tne  oomolete  spatio-temporal  benavicr  of  s^C'-t  )i:nt  p.i-.es  t  **a  v  **»■•.!" ;  optical  -ater'ais  _3j.  It  ..as 
designea  to  be  applicable  ror  any  conceivable  imvii  iser  pea-  cro*if»;  ana  beam  snape  'sc-seg 

into  the  material ) ,  material  distribution  e.i.  joi^on  or  la/e^ej,  aterials  ,-mn  re'r»ctive  'rj-x 
interfaces*  clad  optical  fibers  and  coatings  on  s^ostra'es  . .  "ne  .'.articular  experimental  arrangement 
simulated  in  tne  present  work  is  scne-at ica i 1 j  snown  :n  -iqure  ).  . :  represents  tne  configuration 

employed  py  Smith  et  al.  [9*12]  in  them  studies  of  the  janaqe  tnresnolds  in  alkali  nalioes  at 
\  =  532  nm.  The  FVIHM  pulse  lengtn  is  21  psec. 

i  1  1 

In  the  model  of  the  pulse,  tne  energy  is  carrier  py  10  macroscopic  "photons"  of  10  photons 

each  wnose  trajectories  are  followed  tnroucn  the  medium  1.13].  *  Saussian,  clipped  at  tne  i  ev  .a  lues 
of  the  peak  flux  profile  in  radial  direction  ana  in  time,  is  attributed  to  tne  pulse  ana  botn,  tne 
pulse  and  the  medium,  are  considered  rqtationally  Symmetric.  These  pnotons  ‘  are  injected  v.to  the 
solid  with  an  appropriate  group  velocity  to  mimic  '•efraction  on  tne  solid  interface  ana  *ocusino  b, 
the  external  f  *  2.5*1  cm  lens.  In  tne  absence  of  nonlinear  pnoton-material  interaction,  tne  ‘ocai 
spot  is  located  0.2  cm  behind  the  interface  air— *Jaii.  Aoorooriate  equation  of  motion  Snell's  law  of 
refraction  and  changes  in  group  velocity)  of  eacn  'proton'  are  solved  to  -ma  the  changes  in  direc¬ 
tion,  group  velocity  and  total  proton  content  .correction  for  a&sorotion).  Simultaneously,  tne 
material  rate  equations  <2-5/  are  solved  to  determine  index  n  and  temperature  r.  This  :s  achieved  by 
treaf'ng  all  equations  as  finite  differences  equations  witn  a  time  step  _t.  For  this  purpose  tre 
interaction  volume  is  divided  into  ceils  faming  tne  aria  scnematica 1 1 y  depicted  in  -’jure  2.  're 

cell  size  used  in  the  present  calculation  is  :.r  =  0.9  _n  in  radial  direction  ano  .a  *  7  .m  in  t"e 

direction  of  beam  propagation  along  the  axis,  'he  time  step  cnosen  is  .it  s  0.03489  psec.  A  total  of 

30  cells  in  radial  and  IOC  in  x-direction  form  the  grid.  Further  details  of  "PULSE'  are  given  in 
reference  13. 

Selected  flux  distributions  obtained  from  both  models  are  presented  in  figure  3.  They  show  the 
passage  of  a  21  psec  pulse  througn  tne  ‘'ocal  region  in  \aCl  in  successive  tine  steps  o*  220  it  *$ee 
indication  on  the  right  side  of  each  picture).  The  oea*  rms  field  strength  of  this  pulse  is 
12.7  MV/cm  which  corresponds  to  tne  measured  damage  thresnold  [9,12].  ’ne  numbers  on  tne  left  side 
and  on  the  x-axis  indicate  dimensions  in  -jn.  'he  expected"  locat’on  cr  tne  focal  spot  is  at 
(F,r)  *  (2000  -n,  Q».  'he  direction  of  the  beam  is  from  left  to  rignt.  'he  coior  coce  r  tne  local 
photon  fluxes  is  jiven  in  table  2.  Severe  beam  deformation,  particularly  oast  tne  p  ane  througn  - 
perpendicular  to  the  beam  axis,  is  obvious  for  both  models,  ret,  tne  details  are  distinctly  oi^erent. 
This  feature  is  interesting  because  it  points  to  possibilities  of  discriminating  between  tne  two 
models  of  laser  damage.  For  example,  the  time- integrated  spatial  flux  distribution  in  a  plane  past 
the  focal  volume  could  be  measured  and  compared  witn  calculations. 

The  spatial  distribution  of  tne  carrier  density  at  t  *  1200,  _t  *  41.9  psec  after  arrival  of 

the  pulse  at  (f,r'  *  (1500  „m,  9)  is  snown  in  figures  4  and  5  ana  the  color  code  is  qiven  in  table  5. 

Clearly,  the  spatial  gradient  of  the  refractive  index  associated  witn  sucn  a  distribution  is  respon¬ 

sible  for  the  observed  beam-scattering  away  from  tne  center  of  tne  focal  volume. 

The  spatio-temporal  evolution  of  the  lattice  temperature  T  was  monitored  as  well.  ->"ar*ab)/, 
it  did  not  increase  more  than  4k  above  tne  initial  temperature  of  MC».  a:  any  point  •«  tre  •-•e'-act’jn 
volume  for  either  one  of  trie  two  models.  Thus,  Jamaqe  cel  ievec  to  be  tne  result  :f  m,v:  te~rei'jt-*c 
pile-up,  approaenmg  the  melting  point  Tv  *  1074  is  prevented  in  tne  discussed  case  «m:t 

breakdown  was  observed  in  the  experiment  by  Smith  et  al.  [9,12].  Tms  "ay  only  oe  e^'a^ec  as 
follows: 

i)  Neither  one  of  the  two  models  of  laser  damage  does  account  -or  tne  -eas-’ea  :a-a  :ata  si 
lated  here,  and  quite  possibly,  for  all  other  classic  e*periments  in  j'mI:  naities  *r\.-  a.*  pee' 

used  as  the  basis  for  the  establishment  mainly  of  tne  avalanche  model,  or 

if)  not  the  intrinsic^damage  properties  of  tne  material  nave  been  observe  **  t**e  -*;•  •*••■•  •v 

performed  by  Smith  et  al.  [9,12],  but  sere  extrinsic  •eature  caused  tne  -ate'-iai  */.  *a . - 

intense  photon  field  of  tne  laser  puise. 

Twc  avenues  of  further  -esearen  nave  to  oe  Pursued  tb  resolve  tms  issue.  *  »»  *  '•  *  .-*• 
cerns  the  improvements  of  tne  present  models,  for  example,  pnmar'v  :e*ect  •• 

snown  to  drastically  influence  tne  empties  of  tne  processes  jeeu^r*".-;  :»  i  e  ’ 

In  particular,  the  densities  of  free  carriers  -egutred  ‘or  prea*  town  * s  „e  iM-t  ■»:•••:  * 

order  of  magnitude  smaller  tnan  those  obtained  wnen  tne  snotbcne,*"'CJ  1  r‘M  ic. '  tons ,  -.  .  *■»  •* 

defects,  are  disregarded  as  was  the  sase  '  -  tne  present  mvesti  :a*  ion .  '•**?  :*'♦*  •  "•<• 

ti  on  of  improved  wide  gap  materials  jr.d  tne  :er'*'0'Tiance  or  new  exper’^ents  :n  an  it  tempt  to  - 

mentally  observe  true  intrinsic  breakdown  o  r.e  none  na . 
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Table  1.  List  of  symbols. 


Table  2.  Color  code  for  the  photon  flux 


"c 

:  free  carrier  density,  concentration  of 

[phot,  cm' 

•2sec-'] 

electrons  in  the  conduction  band  [cm*3]. 

<1027 

Black 

n 

V 

22  -3 

1  2.24  *  10  [cm  ]:  concentration  of 
valence  electrons 

!027 

White 

uJ 

:  angular  frequency  of  photons 

1027'5 

Blue 

F 

2 

:  [jhoton  flu*  [number  of  photons  oer  cm  and 

1028 

Green 

sec] 

1028-5 

Red 

A 

*  532  [jn]:  photon  wavelength 

:o29 

Yellow 

E 

:  ms  optical  field  strength  [MV/cm] 

io29-5 

Magenta 

**i 

:  avalanche  ionization  rate  [sec'1] 

1030 

Cyan 

*A)U 


1.45  x  10“"*  [cm3  sec  j  four-photon 
absorption  cross  section 
-17  2 

1.5  *  10  [cm  ]:  polaron  absorption 
cross  section 


Table  i.  Color  code  *'or  carrier  ;ens:t./ 
[percent  of  valence  electron 
density] . 


m 

*  0.75  m  ■  mass  of  free  carriers 

e 

<:!0*5 

.‘lack 

- 

v  2.165  [g  cm’3]: 
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[esu]:  nonlinear  refractive 
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Figure. 3.  Photon  flux  distributions  at  various  times  (in  multiples  of  At 
as  indicated  at  the  rinht  of  each  frame)  in  the  focal  volume.  The  oictures 
in  the  left  column  were  obtained  from  the  avalanche  model,  those  in  the  rirht 
column  from  the  multiphoton-polaron  model.  The  color  code  is  given  in  tabie  2. 


Figure  4.  The  spatial  distribution  of  the 
free  carrier  concentration  reached  at  4.19 
psec  =  1200  x  At  by  avalanche  ionization. 
The  color  code  is  given  in  table  3. 


Figure  5.  The  sDatial  distribution  of 
the  free  carrier  concentration  reached 
at  41.9  psec  by  multiohoton  ionization. 
The  color  code  is  given  in  table  3. 
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IV.  The  contribution  of  spherical  aberrations  to  the 
vestige  structure  induced  by  laser  damage 

Garnett  W.  Bryant  and  Ansgar  Schmid 
Department  of  Physics,  Washington  State  University 
Pullman,  Washington  99164 

Abstract 

The  contribution  made  to  the  propagation  of  a  Gaussian  pulse  by  the 
spherical  aberrations  introduced  by  the  lens  used  to  focus  the  beam  into  the 
test  sample  is  considered  as  a  possible  explanation  of  the  multi-vestige 
structure  observed  in  laser-induced  damage  of  transparent  materials.  The 
paraxial  approximation  for  the  wave  equation  is  used  to  determine  the  inten¬ 
sity  of  the  pulse  on  the  beam  axis  when  the  aberrations  ere  present.  The 
inclusion  of  spherical  aberrations  modifies  the  form  expected  for  a  diffrac¬ 
tion  limited  Gaussian  bean,  shifting  the  main  peak  in  intensity  away  from  the 
focus  and  suppressing  it.  More  importancly,  oscillations  are  introduced  in 
the  intensity  prior  to  the  focus.  Although  the  spatial  arrangement  of  the 
peaks  in  the  intensity  appears  consistent  with  some  of  the  experimental 
results  for  the  vestige  structure,  the  spacing  between  peaks  Joes  not  corres¬ 
pond  to  the  observed  spacing  between  damage  sites.  These  findings  and  oLi.ei 
possible  explanations  of  the  vestige  structure  are  considered  critically. 


PACS  numbers:  42.30.Fk,  42. 60. Kg,  61.50-x,  79.20.Ds 
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I.  Introduction 

One  of  the  intriguing  aspects  of  recent  laser-induced  damage  studies 
on  wide  gap  transparent  materials  is  the  observation  of  distinct  damage 
vestiges  along  the  beam  axis  in  addition  to  the  primary  damage  at  the  focal 

spot.  Such  vestiges  have  been  seen  by  Danileiko^"  at  2.76y,  by  Smith  at 

2  3  4 

1.06u  and  0.532y  and  by  Anthes  and  Bass  also  at  0.532y.  However  the 

observed  damage  morphology  is  not  the  same  for  all  experiments.  Using 

pulses  well  above  the  critical  power  for  self-focusing,  Danileiko  found 

vestiges  separated  by  roughly  50y  and  approximately  10y  in  length.  This 

contrasts  with  Smith's  observations,  using  pulses  below  the  critical  power, 

of  ly  vestiges  separated  by  10y.  Furthermore,  Smith  found  that  the  damage 

sites  were  statistically  spaced  around  the  focal  point.  In  contrast  to  this 

Anthes  found  that  all  the  vestiges  occurred  prior  to  the  focus. 

2  3 

The  avalanche  ionization  model  of  breakdown  has  been  used  ’  to 
explain  the  statistical  occurrence  of  these  damage  sites.  Such  an  explanation 
should  contradict  the  findings  of  Anthes  that  all  the  vestiges  occurred  prior 
to  the  focus.  However  Anthes  also  used  a  pulse  well  above  the  critical  power. 

As  a  result,  his  results  would  be  consistent  with  Smith's  interpretation  pro¬ 
vided  the  strong  self-focusing  present  in  his  experiment  selectively  enhanced 
the  intensity  prior  to  the  focus  and  suppressed  the  formation  of  vestige  sites 
past  the  focus. 

In  contrast,  Danileiko  has  recently  suggested  that  spherical  aberrations 
of  the  Gaussian  pulse  induced  by  the  lens  used  to  tightly  focus  the  beam  into 
the  test  material,  could  significantly  modify  the  intensity  of  the  pulse  near 
the  focal  point,  actually  shifting  the  focal  point  and  adding  a  strong  oscil¬ 
latory  modulation  to  the  intensity  expected  for  a  diffraction  limited,  focused 
Gaussian  beam.  These  strong  oscillations  in  the  intensity  along  the  axis 


2  0 

result  in  regions  of  enhanced  intensity  that  could  cause  the  vestiges.  Pre¬ 
viously,  such  an  explanation  was  used^  to  successfully  explain  the  appearance 
and  characteristics  of  damage  vestiges  in  breakdown  studies  of  gases. 

Danileiko  did  not  give  any  details  of  his  calculations  of  the  effects 
of  spherical  aberrations.  Moreover,  he  restricted  his  attention  to  pulses  at 
2.76p  and  10. 6p.  We  have  performed  detailed  calculations  to  extend  his 
results  to  other  wavelengths  of  interest.  In  addition  we  have  investigated 
the  effects  that  variations  of  the  focal  length  of  the  focusing  lens,  the 
index  of  refraction  of  the  damaged  material,  and  the  position  of  the  lens 
relative  to  the  sample  surface  might  have  on  the  contributions  of  the  aber¬ 
rations  >  We  are  compelled  to  consider  such  variations  to  compare  our  results 
with  experimental  findings  obtained  for  a  variety  of  situations. 

In  Section  II  we  describe  the  calculation  of  the  effect  of  spherical 
aberrations  on  the  propagation  of  a  Gaussian  pulse.  In  Section  III  the 
results  are  presented  and  an  approximate  expression  is  derived  which  predicts 
the  spatial  distribution  of  the  regions  of  enhanced  intensity.  Finally,  in 
SectionlV  we  relate  our  results  to  the  observed  damage  morphology  and  recon¬ 
sider  the  other  explanations  given  for  the  vestiges. 

II.  Theory 

We  use  the  paraxial  approximation  to  obtain  a  tractable  wave  equation 
for  the  field  E(r,t)  which  includes  diffraction  of  the  beam.  Specifically, 
we  assume  that  the  field  can  be  written 

E(r,t)  =  e(z,  xx)  exp(i(kz  -  wt))  .  (1) 

e  is  a  slowly  varying  envelope  function,  z  is  the  position  along  the  beam 
axis  and  x±  is  the  position  perpendicular  to  the  axis.  In  the  paraxial 
approximation^’ ^  we  solve 


0 


(2) 


1  +  2k  Vxe  = 

3z  2k  -t 

with  k  *  um/c  where  n  is  the  local  index  of  refraction.  The  field  can 
be  evaluated  at  any  position  using  known  values  of  the  field  in  one  plane 
(z  ■  Zp)  perpendicular  to  the  beam  axis;2  in  our  case,  either  at  the  lens 
(Zp  ■  z^)  or  at  the  sample  surface  (zp  =  z^) 

-  2tt1(z-z~)~  d2xl  exP(ik^i  “  x1)2/2(z-Zp))e(zp,x|)  .  (3) 

We  consider  the  typical  experimental  situation  with  the  lens  placed 
at  zT  in  free  space  and  the  sample  surface  situated  at  z  (z  >  z  ).  To 
obtain  a  description  of  the  field  inside  the  sample  (z  >  z  )  we  first  use 
Eq.  (3)  to  find  e(zs«xx)» 

e(VXi>  *  2TTi(z°— zL)  1  exP<ik0<Xl  -  Xx>2/2^s  *  *L))e(*L'*I)  (4) 

with  kQ  =  (jj/c  in  free  space.  We  assume  that  E  is  polarized  perpendicular 
to  z.  This  of  course  is  not  strictly  correct.  When  the  initially  plane 
polarized  Gaussian  beam  is  focused  it  obtains  a  small  component  of  polariza¬ 
tion  along  the  beam  axis.  We  ignore  this  small  component.  Since  E  is 
polarized  in  the  plane  of  the  sample  surface,  E  is  continuous  across  the 
boundary.  Consequently,  the  field  inside  the  material  can  be  found  using 
Eq.  (3)  where  the  known  field  is  at  the  surface.  Thus  for  z  >  Zg, 

e(z,xx)  -  27Ti(z-z  ')  }  d2xl  exP<ik(xI  _  xx)2/2(z  -  Zg))e(Zg,xx)  (5) 

s 

m 

with  k  *  nw/c  where  n  is  the  index  of  refraction  of  the  sample.  Strictly, 
the  continuity  of  E  across  the  surface  implies  that 


e(zs  +  0+,  xx)  -  exp(i(kQ  -  k)Zg)£(Zg  -  0+,  x^) 


where  0  is  a  positive  infinitesimal.  The  presence  of  this  constant  phase 
factor  has  no  effect  on  any  of  the  calculations  and  is  ignored. 


Using  Eq.  (3)  to  evaluate  £(z,x^)  we  obtain 

e(z,xx)  *  - - - - -  f  d2x|d^x^  exp(ik(xi-x^)2/2(z-zs)) 

(2iri)  (z-zs)  (zs-zL)  J 

x  exp(iko(x|  -  x^)2/2(zs  -  zL>)e(zL>x^)  . 

The  integral  over  x'  can  be  performed  for  z  >  z  to  yield 

^  r 

e(2**i)  *  2fi(z-~z~)  J  d2xl  exP(ike(^i  "  ^1)2/2(z  -  2L))e(zL,x;) 

with 

_  “M*  -  V 

*  *  Q(ZS  ’  ZL}  +  2  "  ZS  ‘ 


(6) 


(7) 


(8) 


This  expression  for  e  has  Che  same  form  as  Eq.  (3)  except  that  k  is 
replaced  by  k  for  z  >  z  . 

6  j 

Note  that  k  can  be  rewritten  as 
e 


where 


(9) 


z  -  z/n  -  (1  -  n)Zg/n  .  (10) 

The  only  effect  of  the  sample  is  to  change  the  scale  for  z  inside  the  sample. 

The  field  for  z  >  zg  is  the  same  as  that  at  z  in  free  space  provided  that 

z  and  z  are  delated  by  Eq.  (10).  Solving  for  z  we  see  that  z  ■  nz  + 

(1  -  n)Zg.  Thus  the  length  scale  spreads  by  a  factor  n  inside  the  sample. 

We  are  interested  in  the  situation  with  a  focusing  lens  present.  If  Rq  is 

the  focal  length  of  the  lens,  then  the  effective  fr  .1  length  R  when  the 

d> 

sample  is  present  is  R,  +  z.  =  n(R  +  z,  )  +  (1  -  Thus 

&  o 


R  •  R  +  (n-1) (R  +  z.  -  z„).  The  second  term  is  just  the  shift  in  the 
go  o  L  3 

focus  obtained  using  geometrical  optics  to  describe  refraction  at  the 

surface. 

The  spherical  aberrations  introduced  by  the  passage  of  the  beam 
through  the  lens  are  contained  in  our  choice  of  e(z  ,  x  ).  Following  Born 

g 

and  Wolf  we  use 

2  2 

-►  x  x  , 

e(zL,xx)  =»  Eq  exp(-  2^-)  exp(-iko(^g-  +  Bx^)  .  (11) 

o  o 

The  first  factor  in  Eq.  (11)  gives  the  initial  Gaussian  profile,  the  second 

factor  accounts  for  the  distortion  of  the  wave  front  needed  to  focus  the  beam 

a  distance  R  from  the  lens  and  the  last  factor  describes  the  spherical 
o 

aberrations.  For  B  >  0  rays  originating  off  axis  are  focused  more  tightly 
than  those  near  the  axis. 

We  are  interested  primarily  in  the  intensity  of  the  field  on  the  beam 
axis.  In  the  paraxial  approximation  the  intensity  is  given  by2 

I  -  {z  |  e(z,x±)  |  2  +  2ik  (£*(z»x1)V1£;(z,xi)  -  e(z,xi)^ie*(z,x1)) }  .  (12) 

On  axis  the  last  term  makes  no  contribution  so 


I(z)  =  c  z  U 
s 


(12a) 
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where  cg  =  c/n  is  the  speed  of  light  in  the  sample  and  U  =  n  | e (z » x A) j  /Sr 
is  the  energy  density. 

Using  Eqs.  (7)  and  (11)  to  evaluate  I(z)  we  obtain 


1(2)  - 


2  2 
cnk  E 
*  e  o 

32tt3(z-zl)2 


22  ike  2  2  (X1  +  A) 

d  x.d  x.  exp x  (x,  -  xj 


‘r  ~2  ^^(z-z.)  i 

L 
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2  2 

xi  “  x2  4  4 

"  iko(  2R  +  &(X1  "  X2))] 


(13) 


Making  a  change  of  variables  u  *  and  v  ■  -  x^,  two  of  the  inte¬ 

grations  can  be  performed,  yielding 

2  2 

cnk  E  ,oo.a>  ,2.2. 

I(z)  ”  - —  °—  uvdudv  exp(-  ^  -  - — ) 

32tt(z-z,)  **o'o  4a 

L  o 


“  Jo(2 (I^rj  (1  - 


-j~  -  Be(z-zL)  (u2+v2))) 


where  J  is  a  Bessel  function,  R  =  R  k  /k  and  B  =  Bk  /k  .  Letting 
o  e  o  e  o  g  0  6 

u  *  pcosB  and  v  =  psin9,  the  integral  over  6  can  be  done,  reducing  the 
expression  for  I  to  a  single  integral 


cnk  E2  f 
_ e  o 

1Stt(z-z  )  J 


2  2  keP  (Z'ZL) 

Pexp (-p  /4ao)sin(4-(2_-)-  (1  -  —j— 

_  _ L _ 6 _ 

—  (z-2  ) 

(1  (z-zL)Bep  > 

e 


(z-zL)Bgp  )) 


This  integral,  however,  can  not  be  done  analytically.  Consequent] v ,  careful 
numerical  evaluations  of  Eq.  (15)  to  correctly  account  for  the  oscillating 
factor  and  the  pole  were  performed  to  determine  I(z). 

Equation  (15)  can  be  derived  more  directly  when  no  sample  is  present 
by  using  the  diffractive  ray  tracing  method  of  Tappert.;  Furthermore,  Eq.  (15) 
is  analogous  to  but  not  equal  to  the  expression  used  by  Evans  and  Morgan^  to 
study  the  effect  of  spherical  aberrations  on  beam  propagation  in  free  space. 

We  solve  an  approximate  wave  equation  exactly.  By  using  the  Huygen-Fresnel- 
Kirchhoff  theory  of  diffraction,  Evans  and  Morgan  evaluate  approximately  an 
exact  integral.  Thus  it  is  not  clear,  a  priori,  which  method  gives  a  better 
description  of  the  wave  propagation.  We  have  not  yet  made  any  attempt  to 
compare  the  two  methods.  We  have  however  tested  the  paraxial  approximation 
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calculating  corrections  to  I(z)  that  go  beyond  the  paraxial  approximation 
by  using  the  field  given  by  Eq.  (7).  Such  corrections  change  I(z)  by  no 
more  than  a  percent  and  are  negligible  for  our  purposes. 


All  of  the  beam  parameters  except  B  are  well  known.  However,  although 
some  of  the  lenses  used  in  the  laser  damage  experiments  are  described  as  well 
corrected,  no  explicit  values  for  B  are  given.  As  a  result,  we  have  performed 
calculations  for  a  range  of  values  for  B.  To  obtain  some  feel  for  the  magni¬ 
tude  of  B,  note  that  the  ray  trajectory  for  the  field  leaving  the  lens  is 

given  by  taking  the  gradient  of  the  phase  of  the  field.  Using  Eqs.  (1)  and 

2  4 

(11),  we  find  that  the  phase  is  approximately  k^z  -  kQx^/2Ro  -  kQBx^.  Thus 

^  3 

the  ray  direction  is  k  z  -  k  (x./R  +  4Bx,)x,.  The  spread  A  along  the 

oo  ioJ-  i 

beam  axis  between  rays  originating  on  axis  and  rays  originating  a  distance 

a  from  the  axis  at  the  lens  is 
o 


A 


2  2 

4Ba 

o  o 


1  + 


4Ba2R 
o  o 


(16) 


For  Smith's  experiment" at  1.06  u  the  beam  radius  (at  the  1/e  value)  was 

0.1  cm  and  R  =1.27  cm.  Thus  A  is  6p  if  B  =  0.01  cm  ^  and  is  320p 
o 

if  B  *  0.5  cm 

We  can  make  two  crude  estimates  of  possible  values  of  B  for  one  of 

the  lenses  used  by  Smith.  First,  Smith  estimated  the  contribution  made  by 

spherical  aberrations  to  the  radius  of  his  focal  spot  in  the  transverse 

direction  to  be  2.6u.  A  simple  argument  using  geometrical  optics  shows  that 

the  ratio  of  transverse  spread  in  spot  size  to  the  longitudinal  spread  due  to 

2 

refraction  at  the  sample  surface  is  a  /2R  n  .  For  NaCl  with  a  =  0.1  cm 

o  o  o 

and  R  ■  1.27  cm  this  ratio  is  1/60.  If  the  ratio  of  transverse  to  longitu- 
o 

dinal  spread  in  spot  size  due  to  spherical  aberrations  is  in  anyway  comparable 
then  A  ~  160u.  This  implies  a  B  of  0.25  cm  \ 


26 


8  9 

We  can  also  use  expressions  given  by  Born  and  Wolf  and  Ireland  to 
evaluate  B.  If  we  assume  that  the  lens  used  by  Smith  was  a  thin  planoconvex 

_3 

lens  with  a  focal  length  of  1.27  cm  then  B  is  about  0.1  cm  ,  comparable  to 

-3 

the  previous  estimate.  We  perform  the  calculations  for  B  from  0.01  cm  to 
-3 

0.5  cm  .  For  smaller  values  of  B  the  spherical  aberrations  have  little 
effect  except  when  aQ  »  0.1  cm.  Moreover  larger  values  of  B  are  probably 
too  big  to  be  reasonable. 

111.  Results 

_3 

The  results  of  our  calculations  for  X  =  1.06p  and  0.53p,  0.01  cm 
-3 

^  B  ^  0.1  cm  and  0.1  cm  a^  ^  0.3  cm  are  shown  in  Figs.  1-4.  In  each  case 

the  focal  length  Rq  was  1.27  cm  corresponding  to  the  focal  length  of  the 

lens  used  by  Smith.  All  intensities  are  plotted  using  the  same,  arbitrary 

2 

scale  proportional  to  Eq.  In  addition  the  results  are  shown  only  for  propa¬ 
gation  through  free  space.  Results  for  any  other  combination  of  focusing  lens 
and  sample  with  n  ^  1.0  can  be  obtained  from  the  present  results  by  using 
two  scaling  relations.  If  we  write  all  distances  z  in  terms  of  distances 

Az  from  the  geometrical  focus  in  the  sample,  z  =  Az  +  z  +  R  ,  and  corres- 

L  g 

pondlngly  in  free  space,  z  -  Az  +  z  +  R  ,  then  Eq.  (10)  requires  that 

L»  O 

Az  »  Az/n  .  (17) 

Thus  for  a  sample  with  n  i  1.0  the  form  of  the  intensity  on  the  axis  is 
Identical  to  that  for  propagation  through  free  space  provided  that  distances 
from  the  focus  are  rescaled  by  n.  Note  further  that,  provided  the  distances 
Az  are  measured«f rom  the  focal  spot,  the  intensity  distribution  is  indepen¬ 
dent  of  the  separation  between  lens  and  sample. 

Furthermore,  when  n  ”  1.0  the  expression  for  I(z)  is 
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ck  e2  r  exp(-p  /4a0)sin(—  fe-  -  r  • Bp  » 

I(z)  -  - 9-— - -r  [dp  p  - r - - - 2 -  •  (18) 

16tt(z-2l)  Jo  (—  -  —  -  Bp  ) 

L  o 

When  the  results  for  different  focusing  lenses  with  focal  lengths  Rq  and 

are  compared,  the  intensity  at  z*  in  the  system  with  the  second  lens 

has  the  same  form  as  the  intensity  at  z  if  1/ (z-zT )  -  1/R  =  l/(z'-zt)  - 

L  o  L 

1/R'.  For  small  displacements  Az  and  Az ’  from  the  corresponding  focal 

spots  Az*  ~  Az  /Rq  and  the  primary  effects  of  a  change  in  Rq  are  a 

change  in  scale  of  z  and  a  change  in  scale  of  I.  Thus  the  results  for  all 

R*  and  n  can  be  obtained  from  one  set  of  calculations  for  fixed  R  and 
o  o 

n  ■  1.0. 

Several  general  features  appear  when  the  spherical  aberrations  are 
included.  First  the  primary  peak  in  intensity  no  longer  occurs  at  the  geo¬ 
metrical  focal  spot,  but  instead  occurs  prior  to  it.  Moreover,  the  main  peak 
is  smaller  when  the  aberrations  are  included  than  it  is  for  the  corresponding 
Gaussian  beam,  indicating  the  extra  transverse  spread  of  the  beam  when  the 
aberrations  are  included.  This  decrease  can  be  as  much  as  two  orders  of 
magnitude  when  the  aberrations  are  large.  When  the  spherical  aberration  is 
sufficiently  strong  and  the  beam  not  too  narrow,  oscillations  occur  in  the 
on-axis  intensity  resulting  from  constructive  and  destructive  interference. 

We  will  consider  these  oscillations  in  intensity  in  Section  IV  as  a  possible 
explanation  for  the  multiple  damage  sites.  The  intensity  is  only  shown  prior 
to  the  focus  in  the  figures.  Beyond  that  poii.t  there  are  no  more  oscillations 

m 

in  intensity  and  the  intensity  approaches  that  for  the  diffraction  limited 
Gaussian  beam  without  aberration. 

Comparing  Figs.  1-4  shows  the  effect  of  the  beam  radius  aQ,  wavelength 
X  and  aberration  B  on  the  pulse  propagation.  First,  as  the  beam  becomes 


broader,  the  spherical  aberration  (proportional  to  x^)  becomes  more  impor¬ 
tant.  For  small  a^  the  oscillations  disappear  and  the  main  peak  approaches 
in  magnitude  the  intensity  of  the  Gaussian  beam  at  the  focus.  For  large  a^ 
the  main  peak  is  much  smaller  relative  to  the  Gaussian  focal  intensity  while 
the  secondary  peaks  are  much  larger  relative  to  the  main  peak.  The  absolute 
magnitude  of  intensity  increases  for  increasing  aQ  simply  because  a  broader 
beam  is  focused. 

When  A  is  varied  the  aberrations  become  more  important  for  smaller 
A.  This  is  evident  from  Eq.  (18).  Decreasing  A  (increasing  kQ)  is 
equivalent  to  simultaneously  decreasing  Rq,  which  compresses  the  oscillations, 
and  increasing  B,  which  enhances  the  magnitude  of  the  oscillations.  Conse¬ 
quently  similar  trends  occur  for  decreasing  A  that  occur  for  increasing  a^. 
The  main  peak  decreases  relative  to  the  peak  of  the  Gaussian  beam  while  the 
side  peaks  become  more  important.  Although  the  positions  of  the  secondary 
peaks  are  insensitive  to  changes  in  aQ,  they  are  very  sensitive  to  changes 
in  A.  The  pattern  is  compressed  for  small  A.  This  also  occurs  for  a 
Gaussian  pulse  without  aberration.  Pulses  for  smaller  A  are  much  more 
sharply  focused  because  diffraction  is  less  important. 

When  B  is  increased  the  pattern  spreads  out  and  the  side  peaks 
become  more  important.  This  contrasts  with  the  change  that  occurs  when  A 
is  increased  and  the  pattern  spread  out  but  the  side  peaks  become  weaker.  In 
all  cases,  the  trends  observed  for  a^  and  B  at  one  wavelength  occur  at 
all  other  wavelengths  considered  between  0.53y  and  10. 6p.  However  for  very 
long  wavelengths  (A  £  10. 6y)  the  oscillations  are  only  present  for  large  beams 

•  _2 

(a  -  0.3  cm)  or  large  aberrations  (B  ~  0.1  cm  ). 
o 

One  can  obtain  an  approximate  relation  for  the  position  of  the  maxima 
and  minima  in  the  intensity  by  considering  Eq.  (13).  The  position  of  these 


i 


extrema  are  controlled  by  the  zeros  of  the  sine  function  in  Eq.  (18).  The 
zeros  are  given  by  the  solutions  of 


2 

v  -  cv  =  rn7T 


(19) 


with  v  »  k  p2(R  -  z  +  z  )/4R  (z  -  z  ) ,  c  =  (z  -  zt)2R24B/(R  -  z  +  zT)2k 

oo  L  o  L  Loo  Lo 

and  m  an  integer.  The  roots  are 


(20) 


Boots  are  possible  for  all  m  <  ms  such  that  4ircmg  <  1  <  4irc(ms+l). 

When  4irc(ms+l)  *  1  another  zero  of  the  sine  function  occurs.  For  z  slightly 

further  from  the  focus  (z  -  zT  smaller)  c  decreases  and  4rtc(m  +1)  <1.  In 

L  s 

1/2 

the  region  between  the  new  zeros  at  v  =  (1  ±  (1  -  4iTc(m  +1))  )/2c,  the 

s 

sine  takes  on  the  opposite  sign  to  adjacent  regions.  Thus  if  the  sine  is 
negative  between  the  two  new  zeros,  a  new  negative  contribution  is  made  to 
the  integral  of  Eq.  (18)  and  the  intensity  decreases.  Conversely,  if  the  sine 
is  positive  between  the  two  new  zeros,  then  the  integral  will  increase.  Thus 
we  expect  the  extrema  to  occur  for  z  that  solve  4iTcm  =1  (m  >  0) .  These 

HI 

are  given  by 


-  ZL +  <Tsk,1/2'<1  +  ‘Trrr’1'2’  '  <21) 

16?tBR  a 
o 

This  is  not  an  exact  relation  for  the  positions  of  the  extrema  because  it  does 


not  include  the  contribution  made  by  the  other  factors  in  the  integrand. 


Equation  (21)  overestimates  the  distance  between  the  geometrical  focus  and  an 
extremum,  especially  the  distance  of  the  main  peak  from  the  geometrical  focus. 
Equation  (21)  gives  a  much  better  description  ot  the  positions  of  the  secondary 
peaks.  Moreover,  since  the  relative  error  made  for  each  peak  position  is 
similar,  Eq.  (21)  is  much  more  accurate  at  predicting  the  separation  between 


i 


so 


peaks.  Note  that  z^  is  independent  of  a^.  Inspection  of  the  figures  con- 

9 

firms  that  the  peak  positions  are  insensitive  to  beam  width.  Ireland  also 
made  a  similar  observation.  Finally,  Eq.  (21)  confirms  the  numerical  calcu¬ 
lations  that  show  the  peaks  nearer  to  the  focus  to  be  farther  apart. 

IV.  Conclusions 

Although  the  contribution  of  spherical  aberrations  appears  to  be 
Important  for  an  understanding  of  the  damage  morphology  in  gases,  the  results 
for  breakdown  in  bulk  materials  are  inconclusive.  Notice  that  the  spacing 

_3 

between  peaks  when  X  =  1.06u  and  B  *  0.01  cm  (Fig.  1)  is  roughly  20b 

near  the  main  peak.  Moreover  this  separation  must  be  scaled  by  the  index  of 

refraction  (n  “  1.53  for  NaCl).  A  separation  of  30u  is  much  larger  than  that 
2 

observed  by  Smith.  Furthermore,  the  separation  only  gets  larger  for  larger  B 

(Figs.  2  and  3).  Similarly  the  spacing  in  Fig.  4  for  X  *  0.53y  is  much 

larger  than  that  reported  by  Smith.  The  spacing  observed  by  Danileiko  at 

2.76y  is  much  greater  and  thus  more  comparable  to  our  results.  However, 

spherical  aberrations  are  less  important  at  2.76y  than  at  1.06u  or  0.53y  sc 

the  better  agreement  is  puzzling. 

There  are,  however,  several  other  interesting  features  about  the 

observed  morphology.  Anthes  and  Bass  found  no  damage  sites  after  the  focal 

spot.  They  attributed  this  to  a  suppression  of  the  intensity  in  this  region 

due  to  self-focusing  prior  to  the  focal  spot.  However  it  is  equally  possible 

that  no  damage  sites  occur  after  the  focal  spot  because  no  oscillations  in 

intensity  occur  in  this  region.  In  addition,  in  several  of  the  published 

2  3 

pictures  of  damage  morphology  presented  by  Smith,  ’  the  separation  between 
sites  increases  in  the  direction  of  the  beam  propagation.  This  structure 
corresponds  to  our  results  that  the  peaks  are  farther  apart  nearer  the  focus 
provided  all  of  the  sites  observed  by  Smith  occurred  prior  to  the  focus.  It 
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would  be  difficult  to  reconcile  such  a  structure  with  the  random  structure 
predicted  by  the  avalanche  model. 

Another  problem  that  arises  when  spherical  aberrations  are  considered 

is  the  appearance  of  so  many  vestige  sites  for  pulse  powers  just  above 

threshold.  As  the  figures  show,  there  is  a  rapid  decrease  in  the  magnitudes 

of  the  peaks  farther  from  the  focus  except  when  B  or  a^  is  large.  However 

Smith's  results  were  obtained  with  narrow  beams  (a  =  0.11  ±0.01  cm).  Thus 

o 

it  is  not  clear  how  the  secondary  peaks  can  cause  intrinsic  damage  if  the 
main  peak,  is  just  above  the  threshold  intensity. 

Although  the  contribution  of  spherical  aberrations  to  the  damage 
morphology  is  not  yet  compelling,  we  have  not  completed  consideration  of  this 
contribution  when  self-focusing  is  included.  We  expect  the  self-focusing  to 
selectively  enhance  the  peaks  in  intensity,  greatly  magnifying  the  effect 
of  spherical  aberrations.  Even  if  B  were  chosen  sufficiently  small  to  obtain 
the  correct  spacing  between  vestiges,  the  self-focusing  might  enhance  the 
secondary  peaks  sufficiently  to  cause  damage  at  those  sites  as  well  as  at  the 
primary  peak.  We  have  begun  a  study  to  include  the  self-focusing  but  no 
results  are  available  yet. 

Several  questions  still  remain.  First,  even  for  small  aberrations, 
the  main  peak  is  less  than  the  corresponding  peak  of  the  Gaussian  pulse  without 
aberrations.  Such  a  suppression  should  affect  the  interpretation  of  damage 
experiments  but  so  far  has  been  ignored.  Thus,  more  careful  determination  of 
the  spherical  aberrations  present  in  experiments  should  be  made  and  a  more 
accurate  descriptions  of  the  beam  propagation  used.  Even  in  the  simplest 
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case  when,  B  vanishes  the  proper  description  of  the  beam  propagation  should  be 

used.  Smith  used  a  diffraction  limited  Gaussian  beam  with  a  diffraction 

2 

length  k  a  determined  from  the  width  of  the  beam  at  the  focus.  That 


diffraction  length  is  much  shorter  than  the  one  we  use  where  a  is  the 

o 

radius  of  the  beam  at  the  lens.  As  a  result,  the  Gaussian  beam  we  use  is 
much  sharper  near  the  focus  (this  is  primarily  determined  by  how  tight  the 
focus  is  rather  than  by  the  diffraction  length)  than  the  pulse  used  by  Smith. 
Consequently  the  good  success  he  had  in  predicting  the  spatial  range  over 
which  vestiges  could  occur  would  not  be  possible  if  our  diffraction  limited 
Gaussian  pulse  were  used. 

In  summary,  we  find  that  the  contribution  of  spherical  aberrations  to 
be  an  incomplete  explanation  for  the  vestiges.  It  is  unable  with  a  small  B 
to  simultaneously  explain  the  observed  small  spacing  between  vestiges  and  the 
occurrence  of  so  many  damage  sites.  However  the  spherical  aberrations  do 
cause  a  systematic  variation  in  site  distribution  which  is  evident  in  some 
experimental  results.  It  is  still  unclear  how  important  spherical  aberrations 
are  and  whether  self-focusing  will  enhance  their  contribution.  Unfortunately, 
the  explanation  of  Smith  is  also  unsatisfactory.  Use  of  the  overly  broad 
intensity  distribution  allows  the  avalanche  model  to  predict  the  occurrence 
of  damage  spots  over  a  much  wider  spatial  range  than  is  justified.  Much  work 
still  needs  to  be  done  before  a  complete  understanding  of  the  damage  morphology 
is  possible. 
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Figure  Captions 

1.  On  axis  intensity  for  propagation  in  free  space  when  A  =  1.06m  and 

-3 

B  -  0.Q1  cm  .  The  region  shown  is  prior  to  the  focus  z^.  The  solid 

lines  are  for  Gaussian  beams  without  aberration  (a  =0.1  cm  and  0.3  cm 

o 

for  the  lower  and  upper  curves  respectively,  the  maximum  intensity  when 

4 

a  ■  0.3  cm  is  1.7  x  10  ).  The  dashed  curves  include  the  aberrations 

o 

(aQ  *  0.1  cm,  0.2  cm  and  0.3  cm  respectively  for  the  curves  in  order  of 

3 

increasing  magnitude,  the  maximum  intensity  when  a  =  0.2  cm  is  1.3  x  10 

o 

3  3 

and  when  a  =  0.3  cm  is  2.3  x  10  and  1.1  x  10  for  the  first  two  peaks), 
o 

2.  On  axis  intensity  for  propagation  in  free  space  when  A  =  1.06m  and 

-3 

B  *  0.05  cm  .  The  solid  lines  are  for  beams  without  aberration 
(aQ  *  0.1  cm  and  0.3  cm  for  the  lower  and  upper  curves  respectively). 

The  dashed  curves  include  the  aberrations  (aQ  =  0.1  cm,  0.2  cm,  and  0.3  cm 
respectively  for  the  curves  in  order  of  increasing  magnitude) . 

3.  On  axis  intensity  for  propagation  in  free  space  when  A  =  1.06y  and 

-3 

B  =  0.1  cm  .  The  solid  lines  are  for  beams  without  aberration  (a  = 

o 

0.1  cm  and  0.3  cm  for  the  lower  and  upper  curves  respectively).  The 

dashed  curves  include  the  aberrations  (a  =  0.1  cm,  0.12  cm,  0.2  cm,  and 

o 

0.3  cm  respectively  for  the  curves  in  order  of  increasing  magnitude). 

4.  On  axis  intensity  for  propagation  in  free  space  when  A  =  0.53u  and 

-3 

B  *  0.1  cm  .  The  solid  curve  is  for  a  beam  without  aberration 

(a  »  0.1  cm}.  The  dashed  curves  include  the  aberrations  (a  =0.1  cm 

o  o 

and  0.2  cm  respectively  for  the  lower  and  upper  curves,  the  maximum 
intensity  when  aQ  =  0.2  cm  is  600).  The  dotted  curve  includes  aber¬ 
rations  when  a  =  0.3  cm  and  is  scaled  by  10 
o  3 
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Calculations  are  reported  of  the  two-phrtnn  absorption  in  errmamum  using  a  Kane  hand  model.  All 
bands  within  10  eV  of  the  valence-hand  nice  are  considered.  The  mains  elements  are  evaluated  using  li  p 
theory  to  find  the  relevant  wave  functions  The  sniieno.il  n  ' del  of  -Jaidereschi  and  Lipan  is  used  to  ireat 
the  distinction  between  the  Item-hole  and  heasv-nole  h  mds  properlv  Detailed  calculations  arc  made  to  test 
various  approximations  that  simoiity  the  evaluation  i  f  incse  comnlicated  nonlinear  optical  coelficients.  hor 
Ge  it  is  essential  that  the  split-off  hand  he  included  as  an  intermediate  state  and  that  the  nonparaholicny  of 
the  conduction  band  be  treated  correctly  The  results  of  the  most  complete  calculation  are  in  good 
agreement  with  experimental  results. 


I.  INTRODUCTION 

With  the  advent  of  high-intensity  lasers,  the 
study  of  multiphoton  transitions  has  become  in¬ 
creasingly  more  important.  Because  the  selection 
rules  for  single  and  multiphoton  transitions  are 
different,  each  provides  different  information  on  a 
material’s  band  structure.  Multiphoton  transitions 
also  provide  another  avenue  for  studying  exciton 
effects.  Furthermore,  the  study  of  these  proces¬ 
ses  aids  in  the  understanding  and  description  of 
the  propagation  of  intense  electromagnetic  fields 
through  nonlinear  materials. 

In  this  paper  we  present  a  theoretical  study  of 
the  two-photon  absorption  in  germanium.  There 
are  two  important  reasons  for  studying  Ge.  First, 
the  experimental  results  for  Ge  are  not  extensive 
and  those  that  are  available  are  not  consistent.  Z  u- 
bov  et  al.1  measured  the  two-photon  absorption 
coefficient  K,  at  2.36  cm  and  found  it  to  be  1 
cm/MW.  Wenzel  cl  al.1  found  that  K.  was  2.5 
cm/MW  near  2.7  ;un.  This  result  is  inconsistent 
with  the  more  recent  results  of  Gibson  et  tl."  for 
the  wavelength  range  from  2.65  to  3.0  cm.  They 
found  that  Ah  is  less  than  0.75  cm  MW  throughout 
this  region  and  is  onlv  0.50  cm  MW  at  2.7  m. 
Gibson  et  al.  were  one  of  the  first  to  eliminate  me 
contribution  made  to  the  measured  total  aosoro- 
tion  by  those  carriers  generated  in  the  twn-ohnror. 
absorption  process.  Coi'.seoucnrlv.  thev  extracted 
two-photon  absorption  coemcients  that  were  lower 
than  the  previously  measured  total  -losorpti  -n  co¬ 
efficients.  Gibson  et  a!,  also  reanalyzed  earlier 
data1,1  for  InSb  using  their  more  detailed  inter¬ 
pretation  of  the  data.  Their  results  were  umun 
several  orders  of  magnitude  lower  than  the  nre- 
vious  experimental  and  theoretical  estimates. 
However,  more  recent  results'  suggest  that  Gib¬ 
son’s  results  for  InSb  are  too  low  because  me 
conduction  electron  density  it  the  In. so  samples 
was  sufficiently  high  for  electrons  to  fill  the  con¬ 


duction-band  states  needed  for  the  two-photon  ab¬ 
sorption.  We  expect  this  effect  to  be  far  less  im¬ 
portant  in  Ge  with  a  band  gap  five  times  larger 
than  in  InSb.  Consequently  the  Gibson  results 
should  be  the  most  reliable  available  for  Ge. 

The  second  reason  for  considering  Ge  is  that  it 
has  not  been  studied  theoretically  as  much  as  the 
zinc-blende  semiconductors."  ,’1  Only  Arifzhanov 
and  Ivchenko9  have  considered  Ge  in  detail.  They 
provided  a  group-theoretical  analysis  of  the  ex¬ 
pression  for  AT  but  gave  no  explicit  results  lor  A,. 
Ge  has  not  been  extensively  studied  for  several 
reasons.  For  other  semiconductors  one  can  olten 
make  the  approximation  that  the  splitting  A  be¬ 
tween  the  valence  and  split-off  band  edges  is  either 
much  greater  or  smaller  than  the  energy  cap  E. 
between  conduction  and  valence  bands.  Such  ap¬ 
proximations  greatlv  simplify  the  matrix  elements 
that  must  be  calculated.  However,  in  Ge  A  ~ 
Furthermore,  for  the  zinc-blende  semiconductors 
one  can  choose  basis  states  so  that  the  matrix  ele¬ 
ments  involving  states  with  wave  vector  k  trans¬ 
form  simply  when  considering  states  with  difierent 
wave  vectors.  Due  to  the  symmetry  of  the  conduc¬ 
tion  basis  states  in  Ge  ithev  have  the  symmetry 
rr  rather  than  r.‘)  tins  is  not  possible.  These  two 
complications  make  the  evaluation  of  matrix  ele- 
p-.-'-nts  and  the  summation  over  intermediate  states 
more  duficult  to  perform  !  >r  Ge. 

In  Sc-r.  II  we  present  the  calculations  of  e  .  We 
consider  onlv  direct  tw.i-pr.ot.m  t r ins.ti -ms  in  Ge. 
Near  threshold  a  nolo  :  i  created  in  tnc  i.e.iv  -irde 
till. i  or  light-hold  Till  valence  band-,  and  in  elect  r  .n 
fills  a  state  in  the  lowest  conduction  b.uio  <  For 
higher  photon  onerous  we  -ils.i  i •  ■  - 1-. -. i  i ■  ■  r  tr.ir.-iti  ins 
wi, are  a  tvdc  in  the  split-nil  >  o  band  is  cr  ati  i. 
Tv.o-photon  absorption  is  a  two-step  nr-ici ■-■».  'a  e 
thus  consider  all  transitions  where  the  inter:. e- 
iii. .te  state  is  one  n|  tin-  ib  wo  st.it.  -,  or  >  -a  :>■  in 
tile  r  \t  11;  ' he st  conilucti  ill  o.inil  ate.  11. e  ['  hindl. 
T.aese  bands  are  shown  in  Fig.  1.  In  e.icn  case. 
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ENERGY  BANOS 
IN  Ge 


FIG.  1.  Band  structure  of  Ge  including  light-hole 
(Ih),  heavy-hole  (hhi,  split-off  (so),  conduction  <ci,  and 
upper  conduction  (uci  bands.  The  solid  curves  give  the 
effective-mass  approximation.  The  dashed  curves  are 
more  realistic  approximations  for  the  lh  and  c  bands. 


the  bands  mentioned,  we  consider  all  oossible  tran¬ 
sitions  involving  intermediate  states  within  10  eV 
ol  the  valence  edge.  We  use  the  spherical  model 
of  Baldcresehi  and  Lipari* '  to  describe  the  valence 
states.  This  allows  us  to  include  their  degeneracy 
properly  but  not  their  anisotropy.  For  comparison 
we  also  perform  the  calculations  assuming  that  the 
heavy-  and  light-hole  bands  are  degenerate  and  can 
be  described  by  a  single  effective  mass.  This  per¬ 
mits  us  to  determine  the  contribution  of  the  va¬ 
lence-band  structure  to  the  results.  Finally  we 
study  the  effect  of  going  bevond  the  effective-mass 
approximation  for  the  energy  bands.  This  was 
done  recently  for  the  zinc-blende  semiconductors,11 
and  large  increases  in  K',  resulted  when  nonpara¬ 
bolic  bands  were  used. 

In  Sec.  Ill  we  present  the  results.  When  all  of 
the  contributions  are  included  properly  we  obtain 
results  which  are  25'i  lower  then  Gibson’s  results. 
In  view  of  the  uncertainties  in  the  experimental 
results  and  in  some  of  the  parameters  used,  this 
agreement  is  quite  satisfactory.  Exciton  effects 
ignored  in  the  simple-band  model  do  not  appear  to 
make  as  large  a  contribution  as  suggested  bv  Lee 
and  Fan  for  the  zinc-blende  semiconductors.  We 
present  final  conclusions  in  Sec.  IV. 


one  step  is  allowed  (does  not  vanish  for  states  at 
the  center  of  the  band)  and  the  other  is  forbidden 
(vanishes  at  the  center  of  the  band).  There  are  no 
allowed-allowed  transitions  in  Ge.  By  including 


II.  THEORY 


For  light  with  a  frequency  u  and  polarization  fj , 
the  two-photon  transition  rate  is  given  in  second- 
order  perturbation  theory  by 


2»/4t \  V1  f  V'  1*  p  !•/' 1  /  n  *p  i  r>  2  _  _ 

T  n  (ccVmV  J  ,2-)-  Y  t,( lO-t.bO-iiw 


We  do  not  consider  exciton  effccts“,u  in  this  pa¬ 
per  so  the  states  in  En.  > l >  rcter  to  single-particle 
band  states  with  energies  *  sit.  The  sums  are 
over  valence  ami  conduct;  in  states  involved  as 
initial  and  final  states  ar.d  over  intermediate  stales 
J.  Since  we  have  cored  d,n  >ie  transitions,  s  de¬ 
notes  the  wave  vector  a  t-.o  .-fates  <",  ■  .  and 

|7)  Involved  in  the  transit;  in.  /  is  the  oeam  inten¬ 
sity,  t  the  dielectric  function,  and  >n  the  iree-elec- 
tron  mass.  The  tv.-  o-phc.ton  absorption  coefficient 
K is  given  by 


We  assume  that  the  beam  is  random Iv  polarized. 
Thus  our  final  results  are  averaged  over  all  pos¬ 
sible  polarizations. 

We  use  k*n  perturbation  theory  to  describe  the 
wave  funrt:  ts  to  ‘  r  :  ■  -  r  m  T’;:<  is  neces¬ 
sary  because  l..e  ;  i  st.  •:  tt."  p.  >- 


transitions  are  described  bv  the  lirst-order  cor¬ 
rections  to  the  wave  lunctinns.  Tiie  b..sis  states 
for  the  perturbation  the orv  arc  th  >se  at  '•  ■  u.  We 
use  the  oasis  set  ol  Kiweett.'  (lie  matrix  ele¬ 
ments  ot  Eq.  1)  and  those  tieeoeu  in  the  f.  •  p  the  i- 
ry  are  evaluated  using  the  enperimer.t.iih  deter¬ 
mined  matrix  elements  t>  -  O.v  :*  !.  .  ami  »,< 

-0.;-ii22  t Th.e  reil  c  and  so  ti  mds  remain 
doublv  degenerate  lor  i- -- 0.  We  also  treat  th.e  uc 
bands  as  if  they  remain  degenerate  since  their 
contribution  to  K  is  small.  As  a  result,  the  chan¬ 
ges  in  the  t,  so.  and  uc  wave  functions  f  ir  /,• »  0 
can  be  found  using  standard  first-order  perturba¬ 
tion  iheorv.  The  energies  are  given  in  the  cilec- 
tive-mass  approximation  ov  second-order  pertur¬ 
bation  theorv.  To  second  order  the  <•  and  so  t  anas 
remain  degenerate  as  tnev  should  but  the  uc  bands 
•  "dit.  As  nii'iiti  -m  t,  we  i  ii.'i'c  this  - •  ■ . : . »>v 


i  Jill!!*'!  U  S I T '  ’!•-*  ilv  «  I  ^  t  ■  i 


in  i  ...'it'  i 
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TABLE  L  Energv  inns  E-,  and  effective  masses  »>.  nf 
electron  states  in  the  is  c**ucti.  n  icl  and  uncer  cenrhtcrion 
(uc)  bands  and  of  holes  in  'he  nenvv-hol<*  ihhl.  licht-hnie 
(lh),  and  split-off  iso)  bands.  Enemies  are  measured 
relative  to  the  valence-band  maximum  and  are  m  eV. 


a 

SO 

lh 

hh 

c 

UC 

Ea 

0.29 

0.00 

0.00 

0.805 

2.9 

mjm 

0.09 

0.042 

0.33 

0.04 

0.64 

we  list  the  effective  masses  and  energy  gaps. 

The  degeneracy  of  the  Ih  and  hh  states  is  also 
lifted  for  kr  0.  This  cannot  be  ignored  as  it  was 
for  the  uc  band  since  the  valence  states  are  the 
initial  states  of  the  absorption  process.  Conse¬ 
quently,  we  apply  second-order  degenerate  pertur¬ 
bation  theory  to  find  the  correct  combinations  of 
faience  basis  states  that  lift  the  degeneracy  for 
each  C  *0.  When  this  is  done  the  energy  bands  are 
no  longer  spherically  symmetric.  This  makes  the 
evaluation  of  Eq.  (1)  extremely  complicated.  To 
simplify  the  calculation  we  use  the  spherical  model 
of  Baldereschi  and  Lipari10  to  find  the  valence 
states.  In  this  model  the  degeneracy  of  the  valence 
bands  is  accounted  for,  but  the  energies  are  no 
longer  anisotropic. 

In  the  spherical  model  the  effective  Hamiltonian 
of  second-order  degenerate  perturbation  theory  is 

written 

a>_i£>0  (3) 

when  terms  which  lead  to  the  anistropv  are  ig¬ 
nored.  3  is  the  spin  operator  for  a  spin-  ‘  system. 
The  valence  basis  states  if,),  i  -  1  to  4,  that  lift  the 
degeneracy  when,  for  example,  k  =  hi  can  be  writ¬ 
ten  in  terms  of  Fawcett's  states.  In  this  represen¬ 
tation  J,  is  diagonal  and  the  states  < j  =  1.4)  for  •/, 

correspond  to  the  heavv-hole  states  while  the 
states  are  the  light-hole  states.  For  k<*i, 
H  is  diagonalized  bv  a  rotation.  If  •  is  the  angle 
between  k  and  i  and  o  is  liio  azimuthal  angle,  then 
H  is  put  in  diagonal  form  i  I  HE .  by  the  ro¬ 

tation  operators  t  *  =  ex:v/o hi  I  and  U.  -exnt.’O.!,) 
which  diagonalize  the  component  of  J  parallel  to  k. 
The  new  basis  states  o  which  lift  the  degener¬ 
acy  for  wave  vectors  in  the  direction  v  are  given 
by  rotating  the  basis  states  for  h-l\ 

M’jLj  M> 

where >’1,4  for  the  heavv-hole  states  and  is  2  or 
3  for  the  light -hole  states.  Once  t:i.i  .lonrouriatc 
choice  of  basis  suites  1  >r  a  particular  uirection 
of  £  has  been  made,  the  valence  wave  functions, 
including  the  first-order  correction  which  contri¬ 


butes  to  the  momentum  matrix  elements,  are  given 
by  standard  k*p  theory 

!riQ,i,.i)^^±a, 

s  "  £  9 

where  the  sum  is  over  the  c,  uc.  and  so  basis 
states.  In  this  model  the  energy  is  again  given  by 
the  effective-mass  approximation.  The  effective 
masses  are  given  in  Table  I. 

When  we  use  the  effective-mass  energies  the 
kinetic  energy  of  the  electron-nole  pair  created 
by  the  two-photon  absorption  is  «  .  =/;*&"’  2..  . 
where  l/ptv  =  1  Dial  >n.  and  m.  and  »i  are  the  ef¬ 
fective  masses  of  the  conduction  band  and  a  par¬ 
ticular  valence  band.  We  then  get  the  usual  (t  ,)‘  * 
for  the  density  of  electron-hole  states.  Since  the 
amplitude  for  the  two-step  process  is  proportional 
to  fc,  the  transition  rate  from  a  particular  valence 
band  will  depend  on  The  value  of  0  is  fixed  bv 
energy  conservation  and  gives  a  transition  rate 
proportional  to  (uc,)1'-(2'/a)  -  E  \  Because 
(urhJ))s/i  is  four  times  larger  than  (alBl  )'•'-,  the  hh 
band  should  be  much  more  important  in  the  ab¬ 
sorption  process.  Consequently,  it  is  important 
that  we  treat  the  distinction  between  the  two  va¬ 
lence  bands  accurately. 

We  also  consider  the  effect  of  using  realistic 
energy  bands.  The  lower  dashed  curve  of  Fig.  1 
actually  gives  a  better  description  of  the  lh  band, 
as  inspection  of  Fawcett’s  results1-1  shows.  This 
band  does  not  intersect  the  so  band  as  the  effec¬ 
tive-mass  approximation  does  and  the  apparent 
effective  mass  is  larger.  Consequently,  the  use 
of  the  more  realistic  lh  enerev  band  should  in¬ 
crease  the  contribution  of  the  lh  band  to  the  tran¬ 
sition  rate.  A  studv  of  Fawcett's  results  also 
snows  that  the  conduction  band  is  nonparabohe. 
However,  the  modification  he  calculated  does  not 
appear  to  be  as  large  for  the  r  band  as  for  t tie  lh 
band.  Calculations  wall  show,  however,  the'  the 
inclusion  of  the  nonpa rabolintv  '•!  the  conduction 
band  in  the  evaluation  ol  transit:  in  rates  is  nv>ro 
crucial  than  a  proper  ireomrnt  n!  li  e  lh  nand. 

To  obtain  a  »■//#.">/,•  .inuivtic  expression  lor  a 
realistic  lh  band  nr  use  hi  ttieso  calculations,  we 
assume  that  it  has  the  form 


The  factor  was  determined  tiv  requiring 

that  Eq.  (6)  give  the  same  v  correction  as  obtained 
from  fourth-order  perturbation  tt’oorv.  The  other 
constraint  applied  to  d-  term  me  V  and  g ,  wis 
that  the  new  tiaini  have  the  same  >1  ■>;*•  at  hive  /• 
as  the  hh  band.  Tims  1  >/i  -  1  1  «/  ..  Hus 

constraint  is  arbitrary  tint  an  ir.-.iu-ctio:i  ol  the 
bind  i  calculated  bv  F.rvc«  «•  t'  P  .♦  m  n 
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sonable.  We  find  that  AE  =  0.2337  eV  and  Am.„ 

=  0.0481m.  The  results  are  given  as  the  lower 
dashed  curve  of  Fig.  1. 

We  have  not  attemoted  to  determine  an  analytic 
expression  for  a  realistic  c  band  that  would  be 
valid  over  the  entire  range  of  k  in  Fig.  1.  Detailed 
comparisons  of  calculated  K,  with  Gibson’s  results 
are  only  possible  for  tiui  within  0.07  eV  of  the 
threshold.  Moreover,  it  is  not  clear  what  the  ap¬ 
propriate  simple  form  should  be.  However,  in  the 
energy  range  near  threshold  the  following  expres¬ 
sion  is  valid: 


€'(*)=£(5+Ame(feao)')-  (7) 

We  use  the  effective  mass  of  Table  I  for  mc  and 
Fawcett’s  value  Am;  =  -l.l  xlO4  as  an  estimate 
for  the  coefficient  of  the  k'  term.  At  large  k  Eq. 
(7)  is  no  longer  valid  because  it  overestimates  the 
increase  in  the  effective  mass  and  ec  decreases 
for  increasing  k.  Even  for  greater  than  5 
xlO*,  the  bending  of  the  band  shown  in  Fig.  1  is 
probably  too  extreme.  However,  the  comparisons 
we  make  with  the  experimental  values  will  test 
electron  states  for  (k'inY  only  up  to  4.5  xl0\  In 
this  range  we  expect  Eq.  (7)  to  be  an  adequate 
estimate  for  the  mass  enhancement. 

A  comment  should  be  made  about  the  parame¬ 
ters  used.  The  room-temperature  energy  gaps‘J 
are  used  because  the  data  of  Gibson  were  taken  at 
room  temperature.  An  average  energy  is  used 
for  the  uc  band  gap  because  we  ignore  the  splitting 
of  the  uc  band.  The  effective  masses  do  not  ap¬ 
pear  to  be  temperature  dependent.  For  example, 
the  effective  mass  of  the  conduction  band  is  0.04 m 
at  both  1.5  K  (Ref.  15)  and  300  K.1”  Averages  of 
existing  experimental  masses1  are  use<j  for 
m^,,  tn^,  and  muc.  Unfortunately,  the  value  of 
P  that  we  use  is  not  consistent  with  w.  =  0.04w  if 
we  assume  that  the  second-order  perturbation  re¬ 
sult  for  €c  gives  the  correct  effective  mass 


mc  3,‘r  \EC  Ec  -  A  / 


(8) 


This  expression  gives  m.=  0.034m  if  P  =  0.6909 
AVwwa  and  the  room -temperature  energy  gaps  are 
used.  A  value  of  /’  =  0.b24  /■  rn,r-  must  be  used  to 
give  me  «  0.04m  using  Eq.  i8).  Tins  value  cf  P  is 
Significantly  lower  than  tnose  found  cxperirr.entai- 
ly^u.iT  We  perform  calculations  with  this  new 
value  for  P  and  for  m.  =  u.034w  to  test  the  impor¬ 
tance  of  the  uncertainties  in  these  parameters. 

Actually  evaluating  the  sums  in  In.  tl)  is  ex¬ 
tremely  cumbersome  because  a  large  number  of 
terms  are  included  when  the  wave  lunctions  are 
found  to  first  order  m  x-p  tr.eory.  Although  the 
sum  over  conduction  states  in  Eq.  .  li  is  e.ts:i> 


done  because  the  s’"n-up  and  snin-dnwn  states 
made  the  same  c''nrr: 'action,  too  sum  over  valence 
states  and  the  '.ntermevi'.te- state  sum  ore  difficult 
to  perform.  To  proceed,  we  per:  mm  the  inter¬ 
mediate  state  sum  for  fixed  r  and  :.  When  we  in¬ 
tegrate  over  all  directions  of  k  and  average  over 
i),  only  certain  terms  of  the  intermediate-state 
sum  make  contributions.  These  terms  are  evalu¬ 
ated  numerically  and  their  contribution  to  the  tran¬ 
sition  rate  determined.  The  sum  over  r  is  also 
done  numerically.  No  analytic  expression  is  given 
for  T  because  it  would  be  verv  lengthy  and  not 
enlightening.  Instead,  the  final  results  are  given 
in  Figs.  2  and  3.  There  are  several  techniques 
which  make  the  evaluation  of  the  sums  easier. 
These  are  discussed  briefly  in  the  Appendix. 

III.  RESULTS 

All  of  the  results  presented  in  Figs.  2  and  3 
were  found  by  including  the  contributions  of  all 
bands  mentioned  as  intermediate  states  and  by 
making  the  proper  distinction  between  lh  and  hh 
states.  The  solid  curve  of  Fig.  2  was  found  by 
using  the  effective-mass  approximation  for  the 
bands  and  by  ignoring  the  A’  dependence  of  the  en¬ 
ergy  denominators  in  Eq.  (1).  Similar  calculations 
were  performed  to  test  the  importance  of  uc  and 
so  bands  as  intermediate  states.  When  the  first 
was  ignored  K2  decreased  by  ten  to  twenty  percent 


FIG.  2.  Lnergv-disnersion  ettoets  on  two-photon  ab¬ 
sorption.  The  solid  curve  i--  calculated  bv  mciu.i g  all 
bands  but  by  igborin:  toe  t  .l*u-‘nua:.''o  1*1  tne  ryv 

denominators  an  1  nv  u-m.  •  . . :  ■  *  ive-m  a  -  nr  nai¬ 
ra  r..  'I  a<’  w  <rt  o  i  -  and  curve  i  ■  i ■  - 1 ■. vs  V  . . 

d-  :  r  oi  tr.e  denominators,  t:;  ’  I  ,  ,  ■  :  i-  carte  u  —  ^ 

the  realistic  lh  v  md.  i  i  d  ah  0  curve  includes  !>>tn 
efforts.  Tnr  lower  da-  ‘e  d  criv  js  Vo  ih-ir-ri  'n  !rom 
Uc’  so  band  ir-clui  a.;  h  -t.i  ■  't,.. 
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FIG.  3.  Comparison  of  the  experimental  results  of 
Gfbaoo  (Ref.  3)  with  the  simple  band-model  predictions. 
The  solid  curve  is  found  considering  all  effects  except 
the  use  of  the  realistic  c  band.  The  dashed  curve  is 
found  using  the  realistic  c  band.  The  error  bars  are 
tahsn  from  Gibson's  work. 

over  the  energy  range  shown  in  Fig.  2.  However, 
when  both  bands  were  ignored,  h2  dropped  by  a 
factor  of  2.  Ignoring  the  so  band  is  equivalent  to 
assuming  that  A»£c,  hj.  Since  y-\Ec>  this  ap¬ 
proximation  should  not  be  valid  in  Ge.  Although 
ignoring  the  so  band  is  a  good  approximation  for 
some  zinc -blende  semiconductors'  it  is  a  very 
poor  approximation  for  Ge. 

Another  calculation  was  made  to  determine  the 
importance  of  including  the  distinction  between  lh 
and  hh  states.  When  they  are  assumed  to  remain 
degenerate,  the  calculation  of  A",  is  straightforward 
because  the  same  valence  basis  states  can  be  used 
for  all  directions  of  k.  There  is  then  no  need  to 
use  the  spherical  model  to  find  new  basis  states 
for  each  direction  of  k.  When  the  lh  and  hh  bands 
are  treated  as  degenerate  we  can  replace  the  u 
factor  in  the  transition  rate  by  either  v  •  *  •  a) 
or  (V(iicW1  *,-«!),) ;  In  the  first  case  the  results  are 
6%  less  than  the  solid  curve  and  in  the  second 
case,  15'  less.  Thus,  if  we  make  a  fortuitous 
choice  for  an  average  u  v  U.e..  we  do  not  reniace 
it  with  u'M  or  she),  then  ignoring  the  distinction 
between  lh  and  hh  bands  does  not  crucially  aitect 
*r 

In  Fig.  2  we  show  the  effect  of  including  the 
dependence  in  the  cnergv  denominators  (in  the  man¬ 
ner  discussed  in  the  Appendixl  ana  oi  using  trie 
realistic  energy  band  lor  the  light  holes.  As  ex¬ 
pected,  these  contributions  are  only  important  away 
from  the  threshold.  In  the  energy  range  consid¬ 
ered.  the  '•  >rrer:.  .-.sari  mi  inteen  percent.  In¬ 


cluding  the  l!  dependence  of  the  denominators  low¬ 
ers  the  results.  L'sin  t  a  more  realistic  lh  band 
increases  the  contribution  of  the  lh  band  because 
the  lh  density  of  states  ana  the  in  the  matrix  ele¬ 
ment  increase.  The  effects  of  incluaing  the  v  de¬ 
pendence  of  the  denominators  and  using  the  realis¬ 
tic  lh  band  compete  against  eat  h  other.  A  calcula¬ 
tion  whicn  includes  both  differs  little  from  one 
which  includes  neither  in  the  energy  range  up  to 
0.2  eV  above  threshold.  Also  shown  in  Fig.  2  is 
the  absorption  coenicient  for  transitions  starting 
from  the  so  band.  The  shape  is  similar  to  that 
of  the  other  absorption  coefficients.  The  transition 
rate  is  lower  because  the  aensitv  oi  so  states  is 
lower  than  that  of  hh  states.  The  increase  in  en¬ 
ergy  gap  also  decreases  the  transition  rate. 

Different  values  of  P  and  m,  were  used  to  test 
their  effect  on  Ah.  Specifically,  one  calculation 
was  made  with  mc  reduced  to  0.034m  ana  another 
with  P  =  0.624  In  these  calculations  the  t 

dependence  of  the  energy  denominators  was  in¬ 
cluded  and  the  realistic  lh  band  was  used.  The 
new  values  of  P  and  mc  are  too  low  to  be  reliable; 
but,  as  mentioned,  they  are  needed  if  ffq.  (8)  is  to 
give  an  mc  consistent  with  the  P  used.  Since  the 
absorption  coefficient  scales  very  roughly  as  P~ 
and  as  ml/!,  K2  decreases  when  the  new  parame¬ 
ters  are  used.  In  fact,  it  decreases  by  35  '  for  the 
new  value  of  P  and  by  20  '  for  the  new  value  of 
mc.  However,  the  uncertainty  in  K,  should  not  be 
as  large  as  this  suggests  since  the  actual  uncer¬ 
tainty  in  P  and  mc  should  be  much  less. 

A  comparison  between  the  experimental  results 
of  Gibson3  and  the  calculation  that  includes  the 
realistic  lh  band,  the  1’  dependence  of  the  denomi¬ 
nators,  standard  parameters,  and  tile  distinction 
between  lh  and  hh  bands  is  shown  in  Fig.  3.  Our 
comparison  is  limited  because  no  data  are  avail- 
abie  for  >  w  more  than  0.07  eV  from  the  threshold. 
Obviously,  even  when  the  um  ertaintv  in  the  aata 
is  considered,  the  predictions  are  too  imv  bv  a 
factor  of  about  2.  The  agreement  is  idenuate  near 
th.resnold  but  gets  worse  0.o5  oY  hov.r.d  tic  tr.rcs- 
nuld.  The  more  rapid  rise  beyond  O.uij  eV  is  n ■  t 
predicted  by  the  calculation  using  the  simple  . 
band.  However,  when  the  nonpar. ihoi.e  con  ruction 
band  is  used,  a  largo  enhancement  of  the  ••  itouh.ted 
ausorotion  occurs.  In  v.ow  of  the  uncertainty  in 
the  experimental  results.  the  parameters  used  and 
the  approximations  mane,  ti.e  agreement  between 
this  calculation  and  the  experimental  results  is 
quite  satistactorv.  Most  importantly,  the  more 
rapid  rise  in  Ah  at  iagtor  energies  is  evident. 

This  occurs  because  t la  aoo.irent  1  lloetive  mass 
of  the  i  band  increases  nr  increasing  ».•.  itus  ef- 
lict  of  the  nonn.irabolu-.tv  of  tic  «  band  is  es¬ 
pecially  important  nr  tw  >-;  r.  t  n  m  ,n  t:  >:i  -,i 
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scales  roughly  as  the  5  power  of  the  apparent  ef¬ 
fective  conduction  mass  rather  than  as  the  I  power 
as  in  single-photon  absorption.  Consenuently .  the 
two-photon  absorption  should  more  readily  reveal 
the  contribution  of  the  nonnarabolicitv.  This  is  ap¬ 
parent  from  Fig.  3.  Because  we  do  not  determine 
a  realistic  c  band  valid  for  large  k.  we  do  not  de¬ 
termine  the  shape  or  mr  'nitude  of  the  absorption 
coefficient  at  large  energies. 

Zubov1  measured  A",  at  2.36  um  to  be  1.0  em  MW. 
At  this  energy  (0.518  eV)  we  predict  A",  to  be  only 
0.4  cm/MW  when  we  use  the  effective-mass  ap¬ 
proximation.  However,  Zubov  obtained  the  experi¬ 
mental  value  without  correcting  for  the  generated 
carrier  absorption.  As  Gibson  has  shown,  this 
leads  to  A,  values  which  are  too  large.  Moreover, 
use  of  the  nonparabolic  c  band  can  increase  the 
calculated  absorption  coefficient  by  more  than  a 
factor  of  2  for  /?ui  greater  than  0.05  eV  from  thres¬ 
hold.  Thus  the  simple  band  model  should  be  con¬ 
sistent  with  Zubov’s  result  provided  that  realistic 
bands  are  used. 

We  should  note  that  no  attempt  has  been  made  to 
include  higher-order  corrections  to  the  matrix 
elements,  even  though  it  is  apparent  that  higher- 
order  corrections  to  the  effective  masses  are  sig¬ 
nificant.  Unfortunately  too  many  terms  would  have 
to  be  evaluated  to  make  inclusion  of  these  higher- 
order  corrections  feasible. 


are  considered.  For  example,  the  exnton  effects 
increase  the  absorption  in  ZnTe  bv  a  factor  of  10 
and  in  GaAs  by  a  factor  of  2.  In  these  materials 
the  lowest  exciton  nas  a  binding  energy  oi  10  and 
4.4  meV,  respectively,  while  in  Ge  it  has  an  energy 
of  only  1.6  meV.1”  For  this  reason  we  would  fur¬ 
ther  expect  exciton  effects  to  be  less  important 
in  Ge. 

The  effect  on  the  band  states  of  the  intense  fields 
needed  to  measure  two-photon  absorption  has  been 
ignored.  Keldysh.-’'  originally  studied  the  effect  of 
an  intense  field  on  electronic  transitions.  When 
this  theory  was  applied  to  two-photon  transitions.'1 
the  results  were  typically  a  factor  of  10  lower 
than  measured  values  and  were  lower  than  those 
calculated  with  the  band  model  of  Basov.  Thus  it 
would  appear  that  inclusion  of  the  Keldysh  effect 
would  lower  our  results  and  worsen  the  agreement. 
Although  the  simple  band  model  appears  to  be  ade¬ 
quate  for  Ge,  the  importance  of  exciton  and  high- 
field  effects  cannot  be  ruled  out.  These  effects 
require  further  consideration. 
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IV.  CONCLUSIONS 


APPENDIX 

We  must  consider  the  following  operator 


We  have  presented  a  simple  band-model  calcula¬ 
tion  of  the  two-photon  absorption  in  Ge.  The  cal¬ 
culation  treats  the  degeneracy  of  the  valence  bands 
correctly  and  includes  all  allowed  transitions  in¬ 
volving  states  within  10  eV  of  the  valence  edge. 

In  addition,  the  t-  dependence  of  the  energy  denomi¬ 
nators  is  treated  approximately  but  adequately 
and  realistic  energy  bands  for  c  and  lh  states  are 
used.  When  all  of  these  contributions  are  in¬ 
cluded,  the  agreement  between  our  results  and 
Gibson's  results  is  gooa.  The  two  eifects  winch 
must  be  included  to  obtain  reliable  results  are  the 
use  of  the  realistic  <  energy  band  and  tee  use  of 
the  so  band  as  an  intermediate  state.  The  other 
approximations,  ignoring  ttie  v  oenendcnce  of  the 
denominators,  ignoring  the  distinction  between  lh 
and  hh  bands  and  ignoring  ttie  uc  band,  are  all 
more  reasonable. 

Several  effects  have  not  been  considered  in  the 
simple  band  model.  No  attempt  has  oecn  made  to 
include  exciton  eifects  as  done  ov  Lee  and  Km. 

The  good  agreement  tnat  we  get  with  Gibson's  re¬ 
sults  suggests  that  this  meet  should  lie  small. 

Lee  and  Fan  tin  !  a  lar  :e  enh  inrenient  lor  some 
zinc-ble  ■  :e  semic  ;n  luer  irs  w  a  n  exciton  t  lit  rts 


v  .n-p'"Wn-p 

To  determine  T.  we  calculate 

T<c:.W»;r)xr,.W.ic>  • 


(Al) 


(A2) 


We  evaluate  M  for  two  simplified  rases.  In  the 
first  case  we  ignore  the  e  dependence  of  the  ener¬ 
gy  denominators,  using  only  me  energy  runs  to 
evaluate  t  .  v)  -»  ■).  In  >.!us  case  '!  is  the  same 

for  I!i  and  hli  states.  Moreover,  in  the  mt<  rme- 
ui  ite-state  sum.  ttie  sum  over  Hi  and  hh  can  In- 
done  bv  i-'norir.g  the  distinction  between  lii  anti  hh 
states  when  the  non  miiniitnr  is  tne  same  tor  both. 
The  same  valence  basis  states  can  then  be  used 
for  each  wave  vector. 

Wo  also  evaluate  '-/  with  the  c  dependence  of  the 
iieu-  imitator  include-!.  In  this  case  •/  ana  are 
dilferent.  However,  in  the  sum  over  ■  the  sum 
over  lh  and  hh  can  onlv  be  done  sinuilv  il  the  dif¬ 
ference  m  the  two  bands  is  still  len-ired.  L'-ns 
we  in ,ii  e  the  I  iii.iw.i:  •  .ii'nr oramati  -n: 


ma 
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with  r  =  cost-  2)  and  s  =  siniV  2).  The  basis  set  is 
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fer  greatly  from  those  found  when  the  k  depen¬ 
dence  is  ignored  altogether.  Thus  we  expect  Eq. 
(A3)  to  give  an  adequate  approximation  for  the  k 

dependence. 

There  are  several  ways  to  make  the  evaluation 
of  Eq.  (A2)  easier.  First  note  that 

52lklAUi->lJ=  22  52 

*C  M  *€Lt.^  l»€t* 

(A4) 

As  argued  before,  when  we  sum  over  all  lh  and 
hh  states  and  the  operator  is  the  same  for  both,  we 
can  ignore  the  difference  in  lh  and  hh  states  and 
use  the  simpler  valence  basis  states.  Thus  the 
first  term  in  Eq.  (A4)  is  easier  to  evaluate.  In 
the  second  sum  we  must  use  the  lh  states  found 
using  the  basis  states  given  by  Eqs.  (4)  and  (5). 
However,  this  sum  is  easier  to  perform  than  the 
original  one  over  hh  states.  To  see  this  we  need 
to  consider  the  following. 

The  matrix  operators  t'„  and  Ue  take  the  follow¬ 
ing  form”: 


e ‘ 3,4/1  0 
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0 

0 


0 

0 


0 

0 

0 

.-»*»/» 


(A5) 


f/«  = 


c1  -fz^s  vTcsJ  s’ 
-vTc*s  c(l-3.s;)  -s(l-3r')  v3 Vs* 
vTcs2  s(l-3rJ)  c(l-3.s*)  vT<-.< 


V  3CS" 


-  v  3  <•'  s 


(A6) 

where  hh  basis  states  for  k  -  ki  are  jr,)  and  |r.) 
while  jig)  and  |i  .)  are  the  corresponding  lh  states. 
At  other  1c,  the  lh  basis  states  are  given  by  Eq.  pi) 
with  i  =2,  3. 

The  wave  function  is  given  to  first  order  in  k  by 
Eq.  (5).  This  can  be  written  as 


|t',It)  =  s(it)|t',i), 


(A7) 


with 


B  ~f  fl 

The  sum  over  lh  can  now  be  written  as 

52  =  s((t)f  52  k>)k>iVoo. 

V  1=2.3  / 


(A8) 


But 
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Ti  e  product  of  the  rotation  onerators  can  be 
simplified  in  the  following  way: 
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'0  0  O' 

j  0  -k,  -k_/J 2  0 

A=t  (All) 

0  -* **./vT  0 

.0  0  ( l)l/ik ,  0  ^ 

and  ki  =  (kt±i  *,)/ vT- 

The  operator  ,4  defines  new  basis  states  for  the 

Ihband 

for;  =2,3.  (A12) 


These  states  (live  the  same  sum  over  lh  statps 
and  yet  are  simpler  to  work  with.  A  lias  a  simnler 
form  than  b\  .  or  (/_„  and  some  of  the  elements  in 
.1  vanish.  Moreover,  the  elements  of  .1  are  written 
in  terms  of  k.  and  fr,.  When  Eq.  (A2l  is  integrated 
over  all  directions  of  k  only  terms  with  even  pow¬ 
ers  of  kt  and  matching  powers  of  k_  and  k.  will 
contribute.  The  use  of  Eq.  (All)  makes  it  easier 
to  identify  those  terms  in  Eq.  (A2)  that  must  be 
evaluated. 
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I.  Introduction 

Laser  damage  and  laser  annealing  in  semiconductors  have  recently 
lent  a  great  deal  of  motivation  to  fundamental  studies  of  the  inter¬ 
action  of  laser  pulses  with  semiconductors  such  as  Ge,  Si,  GaAs,  etc.  . 

In  addition  to  the  understanding  of  free  carrier  kinetics  and  the  inter¬ 
play  between  heated  carriers  and  the  lattice  background  an  analysis  of 
the  dynamic  laser  self-action  has  now  become  mandatory  if  a  comprehensive 
insight  into  the  intricate  laser  pulse  -  semiconductor  processes  is  to 
be  gained. 

Laser  self-action  can  best  be  described  by  recalling  that  the  dielec¬ 
tric  function  ,  and  in  turn  the  refractive  index,  for  each  material  is 
strongly  dependent  on  the  carrier  distribution.  Thus  a  change  in  this 
distribution  due  to  the  presence  of  an  external  photon  field,  due  to 
changes  in  temperature  or  due  to  other  influences  will  effectively  change 
the  geometrical  optical  properties  of  the  material.  Frequency  shifts, 
self-focusing,  self-ae focusing  and  changes  in  shape  of  laser  pulses  and 
possibly  the  whole  phenomenon  of  laser  damage  in  semiconductors  may  all 
lend  themselves  to  explanations  in  terms  of  such  dynamic  refractive  index 
changes.  We  present  here  a  status  report  on  an  effort  aimed  at  simulating 
such  spatio-temporal  chances  in  Ge  which  is  exposed  to  2.7  Mm  picosecond 
laser  pulses.  We  will  limit  our  discussion  to  temporal  aspects  onlv  and 
will  describe  the  entire  spatial  dynamics  elsewhere.  However,  we  will  be 
able  to  address  and  answer  an  important  question  which  has  recently  been 
raised.  That  is,  Ge  has  been  claimed  to  possess  unusually  high  damage 
thresholds  at  2.7  un  and  10.6  pm.  We  will  comment  on  the  possible  cause 


for  this  effect. 


II.  Damage  Criterium 


From  the  outset  v;e  assume  an  idealized  crystal.  We  disregard  therefore 

any  possible  contribution  to  damage  by  inclusions  or  by  other  defects.  We 

concentrate  instead  on  heating  of  the  phonon  system  which  comprises  an 

Einstein-type  optical  phonon  svstem  and  a  Debve-tvpe  acoustical  phonon 

system.  The  two  systems  are  coupled  to  each  other  and  independent Iv  linked 

to  the  carrier  systems,  i.e.  the  photoexcited  free  carriers  and  holes.  Each 

phonon  system  contributes  to  the  temperature  dependent  mean  square  displace- 

o  2 

ments  of  lattice  constituents,  <u“>  (T)  and  <u  >  (T)  respectively. 

op  ac 

If  these  two  mean  square  displacements  taken  together  become  excessive  (i.e. 
larger  than  the  displacement  sum  at  equilibrium  melting  temperature  T  ) 


<u2>  (T  ) 
ac  ac 


<u2>  (T  ) 
op  op 


> 


<U2>  (T  ) 
ac  m 


<u2>  {TJ 

op  m 


permanent  lattice  modifications  are  expected  and  damage  is  assumed  to  occur. 
Our  purpose  is .therefore , to  perform  calculations  which  permit  monitoring 
of  the  temperatures  in  the  two  phonon  systems  during  the  passage  of  a  laser 
pulse.  Note  that .according  to  our  criterium  the  required  damage  temperature 
has  to  exceed  5600°K  in  the  optical  phonon  system  alone  if  for  some  reason 
the  acoustical  phonon  system  should  nrove  unable  to  significantly  warm  un 
during  the  pulse  Length  of  a  picosecond  pulse.  Note  also  that  the  damaging 
temperature  T  calculated  under  these  circumstances  is  still  an  approxima¬ 
tion  which  will  require  further  refinement.  For  example,  we  presently  eval¬ 
uate  mean  square  displacements  without  accounting  for  anharmonic  effects. 

In  order  to  calculate  the  phonon  system  temperatures  nnd  their  change 
during  the  history  of  a  laser  pulse  detailed  information  about  electron- 
electron,  electron-phonon  and  hole-ntionon  interaction-;  in  -'e  ■;  re. paired. 
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In  the  next  section  we  will  briefly  review  in  form  of  a  block  diagram  how 
the  energy  absorbed  from  the  laser  field  ultimately  ends  un  in  the  optical 
and  acoustical  phonon  systems.  We  will  thereby  specify  some  details  of 
the  different  interactions. 

III.  Energy  Transport 

Photo-excitation  in  Ge  at  2.7u  leads  to  the  formation  of  electron- 

hole  pairs  via  direct  two  photon  interband  absorption.  The  transition  rates 

for  this  process  were  evaluated  in  the  recent  work  by  Bryant  The 

dynamics  of  photoexcited  electrons  and  holes  in  Ge  has  been  studied  in 

(2) 

considerable  detail  by  Elci,  Scully,  Smirl  and  Matter  .  Our  present  model 

uses  their  results  modified  to  take  account  of  the  two  separate  phonon  sys- 
(3) 

terns  and  disregards  Auger  recombination  for  the  time  being.  However,  as 
our  results  will  show,  this  simplification  is  no  longer  justified  since 
appreciable  heating  at  2.7p  requires  very  high  local  free  carrier  concentra¬ 
tions,  a  condition  under  which  Auger  recombination  is  most  effective.  The 
next  phase  of  program  development  will  therefore  include  this  effect  as  well. 

A  further  refinement,  which  is  already  part  of  the  model,  accounts  without 
averaging  for  all  the  conduction  band  valleys  and  allows  for  single-photon 
hole  transitions  between  different  valence  bands  when  the  light  and  heaw 
hole  bands  are  highly  depleted. 

In  our  block  diagram  (Fig.  1)  the  sources  of  energy  gain  from  the  laser 
field  are  listed  on  top.  In  addition  to  two-photon  absorption  we  consider 
free  carrier  absorption  (FG\)  and  free  hole  absorption  (FHA),  which  are  both 
single  photon  absorption  processes  involving  phonons  for  momentum  conservation, 
and  hole  transitions  to  the  split-off  band.  The  last  process  does,  of  course, 
not  generate  any  new  holes  but  increases  their  energv  instead. 


The  respective  incerband  transitions  are  shown  in  Fit;.  2.  Excited 


carriers  in  the  central  conduction  band  valley  are  assumed  to  rapidly 
scatter  into  the  side  valleys  (via  electron-phonon  scattering)  where  they 
thencalize  into  a  quasi-Fermi  distribution  of  temperature  T  .  Note  that 
scattering  into  the  side  valleys  removes  electrons  from  states  which  are 
the  final  states  of  the  two-photon  absorption  process  across  the  band  gap 
and  that  the  build-up  of  carriers  in  the  side  valleys  thus  regulates 
availability  of  those  states  and  therefore  the  transition  rate.  Proper 
accounting  of  the  availability  of  states  demands  as  realistic  a  description 
of  the  bands  as  possible-  Similarly,  the  light  and  heavy  hole  valence 
bands  have  to  be  treated  as  non-degenerate  if  hole  transitions  into  the 
split-off  band  are  to  be  treated  correctly. 

Phonon  emission  and  absorption  leads  to  carrier  cooling  by  intravalley 

relaxation.  Since  the  optical  phonons  have  a  stronger  coupling  to  the 
(3) 

heated  carriers  ,  this  relaxation  will  bring  about  a  slightly  faster 
temperature  rise  in  the  optical  phonon  subsystem  than  in  the  acoustical 
one.  Subsequently  optical  phonons  decay  into  acoustical  phonons  with  a 
characteristic,  temperature  independent  time  constant  .  The  tempera¬ 

ture  rise  in  the  two  systems  is  monitored  in  order  to  detect  damage  according 
to  Section  2. 

The  two  temperature  rises  are  given  by 
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where  optical  and  acoustical  phonons  have  each  been  assigned  half  the  con¬ 
tribution  to  the  total  specific  heat  cv  and  where  x  is  the  optical  nhonon 
decay  lifetime,  NQp(T)  the  optical  phonon  number  at  temperature  T  and  of 
uniform  energy  ,  and  du/dt  the  respective  rate  of  change  in  energy 

density  due  to  energy  transfer  during  free  carrier  absorption  and  intravalley 
relaxation.  These  energy  transfer  rates  depend  on  the  carrier  and  hole 
density  and  respective  energy  distributions  and  therefore  on  the  temperature 
of  the  electron  and  hole  system.  We  will  refrain  from  presenting  the  per¬ 
tinent  expressions  here  since  they  are  rather  involved  and  will  be  published 
when  we  have  arrived  at  a  comprehensive  picture  of  the  laser-semiconductor 
interaction.  It  suffice  to  say  that  the  numerical  procedure  of  modeling  the 
temporal  aspects  of  this  interaction  involves  solving  a  set  of  five  coupled 
differential  equations  each  of  which  contains  integrations  over  the  carrier 
distributions  in  each  valley  and  each  valence  band. 

IV.  Results 

We  have  tested  the  damaging  properties  of  several  2.7jj  laser  pulses 

in  the  picosecond  regime.  Pulse  lengths  of  5,  10,  and  20  psec  were  chosen 
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at  peak  fluxes  of  10  ,  10  and  10  photons  per  cm  and  sec,  respectively. 

For  reasons  of  operational  simplicity  the  Gaussian  temporal  waveform  of  each 

3 

pulse  was  truncated  at  the  1/e  peak  intensity  points.  Initial  conditions 
for  the  carrier  chemical  potential  and  carrier  temperature  were  determined 
by  long  hand  and  the  dependence  of  the  final  results  on  the  precise  choice 
of  these  initial  conditions  was  tested. 

We  present  in  Fig.  3  and  4  the  free  electron  density  N  and  free 

28 

carrier  energy  density  £  per  unit  volume  as  both  evolve  during  a  10 
2 

photons/cm  sec  neak  flux,  20  nsec  pulse.  Without  an  effective  drain  on  the 
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carriers,  such  as  Auger  recombination,  one  finds,  as  expected,  a  noid  in¬ 
crease  in  the  carrier  density  in  the  early  part  of  tne  pulse  without  anv 
significant  reduction  in  the  later  part.  A  minute  reduction  which  is  noti- 
cable  only  numerically  and  not  on  the  graph  of  Fig.  3  accounts  for  two- 
photon  stimulated  emission  which  is  included  in  our  two-photon  transition 
probability.  The  energy  density  E,  of  the  free  carriers  and  holes  shown 
in  Fig.  4  exhibits  similar  behaviour  in  the  first  half  of  the  pulse.  Along 
the  trailing  edge  of  the  pulse  the  energy  density  diminishes  first  on 
account  of  the  two  phonon  systems  and  stabilizes  seemingly  towards  the  end 
of  the  pulse.  This  effect  is  a  model  artefact  and  can  be  explained  in  con¬ 
junction  with  Fig.  5.  There  the  systems'  temperatures  are  plotted  for  the 
free  carrier  and  optical  phonon  systems  as  a  function  of  time.  The  free 
carriers  heat  in  accordance  with  the  instantaneous  availability  of  photons 
and  of  final  states  for  transitions.  The  electron  temperature  follows 
therefore  a  modified  Gaussian  as  one  would  expect  for  a  temporal  Gaussian 
Intensity  distribution  in  the  laser  pulse.  However,  since  the  temperature 
rise  is  accompanied  by  additional  free  carrier  generation  while  the  electron 
temperature  decline  is  lacking  a  corresponding  recombination  process  across 
the  band  gap(the  available  conduction  band  states  rapidly  fill  up  during 
carrier  cooling,  in  the  later  part  of  the  pulse,  thus  keeping  energy  in  the 
electron  system  which  otherwise  could  be  transfered  in  electron-phonon  pro¬ 
cesses.  This  also  explains  why  the  optical  phonon  temperature  rise  proceeds 

at  a  less  rapid  rate  in  the  late  part  of  the  pulse  than  earlier. 

23  2 

We  also  note  that  a  10“  phot/cm  sec  peak  flux,  20  psec  pulse  raises 

the  optical  phonon  temperature  merely  to  552  °K.  Contrary  to  some  earlier 
(3) 

expectations  however,  the  acoustical  phonon  temperature  rise  does  not 
lag  behind  the  optical  one  and  the  two  svstems'  temperatures  agree  to  within 


less  Chan  1%  over  the  whole  pulse  length.  The  optical  phonon  temperature 
does  therefore  not  have  to  rise  to  5600  °K  in  order  for  the  damaping  phase 

transition  to  take  place. 

It  is  tempting  to  compare  these  results  with  earlier  work  done  at 

(2) 

1.06  um  .  While  we  find  qualitatively  plausible  agreement  comparison  of 

numbers  becomes  difficult  in  light  of  the  different  experimental  conditions 

chosen,  in  each  model.  The  work  of  Reference  2  rests  on  the  assumption  of 

a  square  pulse,  for  instance,  while  the  present  work  assumes  a  Gaussian 

temporal  profile.  Since  the  dynamic  processes  are  nonlinear  no  simple 

scaling  relation  between  the  pulse  shapes  can  be  established.  Our  model 

calculations  do,  however,  confirm  that  ideal  Ge  requires  very  high  fluxes 
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even  at  2.7  pm  for  damage  to  occur.  The  flux  regime  of  10  photons/cm  sec 
at  which  we  predict  damage  to  occur  borders  on  fluxes  which  in  much  wider 
band  gap  materials,  at  larger  photon  energies  and  under  physical  processes 
quite  different  from  the  ones  considered  here  also  produce  damage.  These 
processes  have  been  identified  to  be  highly  nonlinear  and  are  thus  requiring 
very  high  fluxes  to  come  about.  In  contrast,  free  carrier  generation  at 
2.7  pm  in  Ge  is  only  a  two-photon  process  and  damage  due  to  these  carriers 
was  intuitively  expected  to  ensue  at  much  lower  fluxes.  Our  calculations 
show  that  for  ideal  Ge  such  expectations  are  unjustified. 

In  the  next  phase  of  our  calculations  several  additional  physical 
phenomena  will  be  included  into  the  model.  Among  them  are  Auger  recombina¬ 
tion  of  free  carriers  and  band  gap  narrowing  effects  as  well  as  the 
dynamic  dielectric  constant  which  will  allow  a  comprehensive  spatio-temporal 
modeling  of  beam  propagation. 

This  work  was  supported  by  AFOSR  through  Contract  -IV  49620-78-C-0095 . 
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Figure  Captions 

1.  Block  diagram  illustrating  the  energy  transfer  in  highly  photo- 
excited  Ge  from  the  2.7  pm  photon  field  to  the  damage  event. 

2.  Band  structure  in  germanium  representing  the  excitation  processes. 

_3 

3.  Electron  density  Ng  cm  as  a  function  of  time  normalized  to  the 

28  2 

pulse  length  t  for  a  peak  flux  of  10  photons/cm"  sec. 

4.  Variation  of  free  carrier  and  hole  energy  density  £  with  time 

2  8 

normalized  to  pulse  length  t  for  a  peak  flux  of  10  photons/ 

2 

cm  sec. 

5.  Variation  of  electron  system  temperature  (left  vertical  scale) 

and  optical  phonon  temperature  T  (right  vertical  scale)  with 

op 

time  normalized  to  pulse  length  t  for  same  pulse  as  in  Fig.  4. 
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ATTACHMENT  A 

Computer  Code  "PULSE" 

PULSE  is  a  two-dimensional  photon-particle  simulation  code  based 
on  RAMSES.^  ^  We  have  adapted  and  modified  RAMSES  to  make  it  suitable  for 
the  calculation  of  the  complete  spatio-temporal  photon-flux  distribution 
as  well  as  that  of  the  photon-induced  material  modifications  which  are 
encountered  when  a  short  intense  laser  pulse  propagates  through  a  nomi¬ 
nally  transparent  medium.  PULSE  is  designed  to  simulate  the  propagation 
for  any  conceivable  initial  laser  intensity  profile,  beam  shape  (e.g.  focused, 
parallel  or  defocused) ,  material  distribution  (e.g.  uniform  or  layered  such 
as  refractive  index  interfaces,  clad  optical  fibers  and  coatings  on  sub¬ 
strates,  etc.)  and  physical  mechanism  of  photon-induced  material  modifica¬ 
tions  (e.g.  photon-induced  changes  in  the  dielectric  response  function,  etc.) 

To  test  the  developed  code  we  have  so  far  concentrated  on  the  simu¬ 
lation  of  an  intense  frequency-doubled  21  psec  pulse  from  a  Nd-glass  laser. 
This  pulse  is  focused  through  a  plane  interface  into  single  crystalline 
NaCl  with  the  aid  of  a  plano-convex  lens  of  1"  focal  length.  The  location 
of  the  geometrical  focal  point  is  chosen  to  be  0.2  cm  behind  the  interface 
inside  the  crystal  (Fig.  1).  This  configuration  is  the  one  used  by  Smith 

(  p ) 

et  al.  * ' in  their  studies  of  the  damage  threshold  in  alkali  halides  at  that 
wavelength. 

In  the  following  we  will  use  this  particular  example  to  describe 
the  features  of  PULSE. 

In  the  model  of  the  pulse,  energy  is  carried  by  a  finite  number  of 
macroscopic  photons  whose  trajectories  through  the  medium  are  the  object  of 
the  calculations  to  date.  A  clipped  gaussian  profile  is  attributed  to  the 
pulse  and  both  the  pulse  and  the  medium  are  considered  rotat iona ! I v  symmetric 
so  that  energy  varies  in  two  dimensions  only.  A  variation  in  the  direction 


of  propagation  corresponds  to  a  time  variation  and  a  variation  perpendicular 
to  propagation  describes  the  radial  properties  of  the  pulse. 
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In  these  present  calculations  the  interaction  between  the  photon 

pulse  and  the  medium  involves  multiphoton  generation  of  free  carriers,  n  , 

which  la  turn  effects  the  refractive  index  n.  In  addition,  a  self-focusing 
2 

term  n^E  is  added  to  n,  where  E  is  the  rms  optical  field  strength.  Local 
variations  of  n^,  n,  and  E  are  calculated  as  functions  of  time  and  the  spa¬ 
tial  coordinate  r  .  The  aim  of  this  simulation  is  to  study  beam 

deformation  of  a  focused  beam  having  a  clipped  Gaussian  flux  profile. 

Initially  the  photons  are  distributed  uniformly  within  the  pulse 
length  and  beam  width.  They  are  directed  to  the  focal  point  and  then 
repositioned  to  take  into  account  refraction  at  the  plane  surface  of  the 

interface  between  vacuum  and  NaCl  (n  =  1.53).  The  clipped  Gaussian 

o 

profile  is  represented  as  follows:  the  intensity  per  cell  "k,j"  given  by 


Here  c  is  vacuum  velocity  of  light,  E  the  total  energy  content,  h^,  h^  the 
cell  size,  R  the  radial  scale  factor  and  2x  the  pulse  length,  a  is  a  free 
parameter  used  for  clipping  the  Gaussian  (1.5  used  to  date). 

For  NaCl  the  local  refractive  index  is  obtained  from  our  previous 
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work 


assuming  optical  phonon  scattering: 
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(b)  0.1  x  0.0022  cm2;  RUN  3 
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F(r)  is  flux  in  cm  sec 

6  -  1.737  x  10-33  (ftui) 

rx^r)  *  (0.5  x  (e^r)  +  f(r)]}5* 

f(?)  =  {ej(r)  +  e2(r) }h 

e.(r)  =  n2  -  s  g(r) 

1  o  w 

C2(J)  '  e20V(;> 

e2Q  ■  5* 

hu  ■  2.33  eV  (photon  energy), 

hoj^  *  .03  eV  (phonon  energy), 

22  -3 

Nvq  *  2.24  x  10  cm  (density  of  valence  electrons) 

s  «  N  x  1.839  x  10_21/(hu)2 
v  vo 


and  g(r)  is  the  ionization  fraction  n  /N 

c  vo 
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Gradients  in  nc(r,t)  cause  gradients  in  n(r,t)  which  lead  to  refrac¬ 
tion.  The  photons  change  directions  according  to 
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where  are  vectors  normal  to  r,  the  particle  directions. 

Since  the  local  velocity  of  light  v  in  the 

g 

the  refractive  index, the  group  velocity  straggling  factor 


medium  depends  on 
0  can  be  computed: 


B(r)  -  v  (r)/c 


n(r)  + 


-  e,(r) 


2n1(r) 


1  + 


e1(r)  -  4  e2Qe2(r) 


f(r) 


The  multiphoton  ionization  ra  te  is  given  by 

R(r)  -  o4(I(r)/fia))4 

-  (saf(;»4 

(S  -  1.951  *  10  for  F  in  previous  units).  Then  the  ionized  fraction 

A 

g  Is  computed  from 

4  -  (l  -  g)R  2) 

so  that  nc  *  8  *  Nv0  • 

The  photon  motion  is  numerically  described  on  a  rectangular  grid 
either  traveling  with  the  pulse  as  in  RAMSES  or  fixed  around  the  focal  volume. 
The  latter  arrangement  was  used  to  produce  the  movie  described  in  Section  IV. 

A  step  AS  is  specified  in  the  input  and  is  ~  1/3  h^,  the  cell  size 
in  x,  so  that  photons  do  not  move  more  than  a  cell  length  in  a  time  step  At, 
where 

At  *  AS/c  .  an  arbitrary  definition. 

Details  of  the  spatio-temporal  resolution  (time  steps,  cell  size  of  the  grid) 
are  given  in  Section  IV  as  well. 

Photon  motion  is  based  on  a  simple  predictor-corrector  scheme.  First  we 
calculate  the  change  in  particle  position  due  to  the  gradient  in  refractive 
index: 
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A^t  .  h 

x  l  y(.3xj  x  [3yJ  J  > 

y  L  y  l3xJ  x[3yj  J  > 


The  variable,  1,  denotes  the  value  of  the  variable  at  time  step  1.  The 
new  directions  at  time  step  (£+1)  are  found  to  be 


T1+1  -  (t1  +  At+\  )-D 
X  X  X 


t1+1  -  Ct8,  +  a*+1t  )*d 
7  7  7 


» ,  mol  ♦  *l4V  ♦  pJ  ♦  iw\)2J 


We  now  include  the  velocity  straggling  factor.  Ar  =  v^  At  so  that 
the  updated  particle  positions  become: 


l  l+l 

1+1  i  ,  n,±^  TX  .  TX  A  f 

x  ■  x  +  $(r)  -y  +  — 2 —  AS 


l  1+1 

-1+1-/  +  ew  f +  V  is 


The  power  density  in  each  cell  is  given  by 


J..  -  cE  W  /V., 
jk  jk  jk 


All  our  work  to  date  consists  of  attempting  to  simulate  the  break-down 

morphology  (which  we  expect  to  be  caused  by  beam  deformation)  observed  by 
(2  ) 

Smith  et  al.  At  the  damage  threshold  a  small  number  (<  4)  micro-sites 

are  found  in  NaCl.  This  goal  dictated  the  particular  choice  of  the  experi¬ 
mental  arrangement  shown  in  Fig.  1. 

In  the  early  calculation  it  was  established  that  a  moving  grid  over 
a  pulse  is  not  practical.  In  order  to  improve  the  spatio-temporal  resolution 
the  code  was  modified  to  use  a  fixed  rectangular  grid.  To  save  memory  space, 
a  fixed  double-triangular  grid  over  the  focal  volume  is  presently  being  experi¬ 
mented  with.  The  resolution  possible  with  this  grid  is  expected  to  be  close 
to  resolving  spatial  features  in  the  order  of  a  few  microns. 

Dnlike  RAMSES  we  now  have  separate  routines  for  the  group  velocity 
and  refractive  index.  Therefore  the  shifting  routine  in  RAMSES  is  eliminated 
and  the  code  is  now  no  longer  restricted  to  gases  but  is  applicable  to  any 
non-linear  (or  linear)  pulse  propagation  problem  where  the  interaction  area 
is  smaller  than  or  equal  to  the  r  -  t  volume  of  the  pulse.  It  is,  further, 
a  simple  matter  to  turn  "on"  or  "off"  the  interaction  between  matter  and 
photons,  which  is  equivalent  to  simulating  classical  (linear)  optics. 
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